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ABSTRACT 
We investigate the origin of the intrinsic scatter in the correlation between black hole mass (Mbh) and bulge 
luminosity (L\, u {) in a sample of 45 massive, local (z < 0.35) type 1 active galactic nuclei (AGNs). We derive Mbh 
from published optical spectra assuming a spherical broad-line region, and Lbui from detailed two-dimensional 
decomposition of archival optical Hubble Space Telescope images. AGNs follow the Mbh -Lbu\ relation of inactive 
galaxies, but the zero point is shifted by an average of AlogMfiH ~ -0.3 dex. We show that the magnitude of the 
zero point offset, which is responsible for the intrinsic scatter in the M^H-Lbui relation, is correlated with several 
AGN and host galaxy properties, all of which are ultimately related to, or directly impact, the BH mass accretion 
rate. At a given bulge luminosity, sources with higher Eddington ratios have lower Mbh- The zero point offset 
can be explained by a change in the normalization of the virial product used to estimate Mbh, in conjunction with 
modest BH growth (^\Q)%-AQ%) during the AGN phase. Galaxy mergers and tidal interactions appear to play an 
important role in regulating AGN fueling in low-redshift AGNs. 

Subject headings: galaxies: active — galaxies: bulges — galaxies: fundamental parameters — galaxies: 
photometry — quasars: general 



1. INTRODUCTION 

Early-type galaxies commonly contain a central black hole 
(BH) whose mass strongly correlates with the bulge luminos- 
ity (Kormendy & Richstone 1995; Magorrian et al. 1998) and 
stellar velocity dispersion (Gebhardt et al. 2000; Ferrarese & 
Merritt 2000). Lower-mass BHs found in late-type spirals and 
spheroidal galaxies follow a similar Mbh _ c* relation (Barth 
et al. 2005; Greene & Ho 2006a) but apparently a different 
MBH-^bui relation (Greene et al. 2008). The BH-host scaling 
relations suggest that BHs play an important role in galaxy for- 
mation and evolution (e.g., Granato et al. 2004; Di Matteo et 
al. 2005; Robertson et al. 2006). Understanding the mecha- 
nism by which BHs coevolve with their hosts impacts current 
models of cosmological structure formation (e.g., Bower et al. 
2006; Croton et al. 2006). 

A key, unanswered question is how the BH-host scaling re- 
lations originated. This issue can be addressed by extending 
the BH-host scaling relations to active galaxies — wherein, the 
BH, by selection, is currently still growing — and by tracking 
the scaling relations to higher redshift to see when and possi- 
bly how they were established. BH masses in type 1 (broad- 
line, unobscured) active galactic nuclei (AGNs) now can be 
routinely estimated to reasonable accuracy (~ 0.3-0.5 dex), 
from the "virial method" using single-epoch ultraviolet or op- 
tical spectra (e.g., Kaspi et al. 2000; McLure & Dunlop 2001; 
Vestergaard 2002; Greene & Ho 2005). More challenging to 
obtain are reliable measurements of the underlying host galaxy, 
particularly of the bulge component, which is maximally af- 



fected by the bright AGN core (e.g., McLure et al. 1999; Floyd 
et al. 2004; Nelson et al. 2004; Greene & Ho 2006b; Kim et al. 
2007), although substantial progress has been made. 

Recent studies present tantalizing evidence that the BH-host 
scaling relations for active galaxies evolve with redshift, even 
by z ~ 0.4 (Woo et al. 2006; Treu et al. 2007) and as far back 
as z ~ 4 (Peng et al. 2006a, 2006b; Shields et al. 2006; Ho 
2007). Compared to local, inactive systems, the sense of the 
evolution is that for a given host galaxy mass or gravitational 
potential, higher redshift AGNs have a larger BH mass than 
similar systems at low redshift. Taken at face value, this sug- 
gests that the growth of the BH precedes, or at least outpaces, 
the growth of the galaxy at higher redshifts. On the other hand, 
local (z ~ 0) AGNs seem to behave quite differently. McLure 
& Dunlop (2002) studied a sample of 72 nearby AGNs and find 
that they roughly follow the same Mbh _ ^bui relation defined 
by inactive galaxies, albeit with a somewhat greater scatter. In 
the Mbh _ c* relation of AGNs, highly accreting AGNs seem 
to have a different normalization in the sense that they tend to 
have a lower Mbh for a given cr* (Greene & Ho 2006b; Shen 
et al. 2008). A similar trend is seen by Ho et al. (2008), us- 
ing H I line widths to constrain the gravitational potential of the 
underlying host galaxy. 

As a concrete step toward establishing a robust z = 
baseline for comparison with high-z studies, this series of 
papers (Kim et al. 2007, 2008; M. Kim et al., in 
preparation) focuses on quantifying the local Mbh _ ^bui 
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FIG. 1. — Distribution of A/bh vs. Eddington ratio. Radio-loud and radio- 
quiet objects are denoted by rilled and open symbols, respectively. 

relation for active galaxies. Among nearby inactive galaxies, 
BH mass correlates only marginally less tightly with bulge lu- 
minosity or mass than with bulge stellar velocity dispersion 
(Marconi & Hunt 2003; Hiiring & Rix 2004; Novak et al. 
2006). Moreover, the M BH - A>ui relation for active galax- 
ies shows no large systematic differences from that of inac- 
tive galaxies (McLure & Dunlop 2002). This suggests that the 
Mbh - ^bui relation can be used as a useful substitute for the 
Mbh _ cr* relation, being an especially effective observational 
tool to track the cosmological evolution of the BH-galaxy con- 
nection (Peng et al. 2006a, 2006b; Treu et al. 2007). Whereas 
stellar velocity dispersions are difficult, if not impossible, to 
measure for distant quasars, for example, photometric mea- 
surements of quasar hosts continue to be feasible even out to 
high redshifts, either through direct imaging (e.g., Kukula et al. 
2001; Ridgway et al. 2001) or through strong lensing magni- 
fication (Peng et al. 2006b). In a companion paper, Kim et al. 
(2008) demonstrate that the bulge luminosity of type 1 AGNs 
can be measured to a reasonable accuracy (~ 0.5 mag) in Hub- 
ble Space Telescope (HST) images, even in the regime when the 
active nucleus far outshines the galaxy. 

Instead of characterizing the full Mbh _ ^bui relation for 
AGNs, this paper restricts itself to only one important aspect: 
the origin of the intrinsic scatter. By choosing a sample for 
which we can estimate reliable BH masses and for which we 
can derive robust measurements of bulge luminosity from HST 
images, our objective is to quantify the true intrinsic scatter of 
the relation and to characterize possible variations of the scat- 
ter with physical properties of the AGN, host galaxy, or envi- 
ronment. By elucidating the physical drivers that influence the 
scatter of the MBH-^bui relation, we hope to gain insights on 
how the BH-host galaxy relations were established. 

This paper is structured as follows. We describe the sample 
selection in §2. We present the image-fitting procedure for mea- 
suring bulge luminosities and our image decomposition results 
in §3. We investigate the Mbh-L^u] relation for type 1 AGNs 
in §4. Finally, Section 5 discusses the origin of the intrinsic 



FIG. 2. — Distribution of redshifts for our sample. The open histogram shows 
the total sample; the hatched histogram shows the objects with low Eddington 
ratio (ibol/^Edd5~ 0.1); the dashed histogram shows the radio-loud objects. 

scatter in the M BH - L^ui relation, ending with a summary in 
§6. Throughout we adopt the following cosmological parame- 
ters: H = 100/t = 71km s" 1 Mpc" 1 , Q m = 0.27, and Vt A = 0.75 
(Spergel et al. 2003). 

2. SAMPLE SELECTION 

Our initial selection begins with all AGNs known to pos- 
sess broad emission lines (type 1 objects), and have reasonably 
deep, and non-highly saturated optical images in the HST pub- 
lic archive. Since we are interested in establishing the local 
MBH-^bui relation, we only consider sources with z < 0.35. 
Next, we carefully search the literature for published measure- 
ments of broad emission-line widths (either Ha or H/3), which 
are needed for calculating Mbh- Although this step is neces- 
sarily somewhat subjective, we try to be consistent in select- 
ing only objects that have line widths with published error bars 
<T0%. The availability of spectrophotometric measurements 
is not essential for us because our BH masses ultimately make 
use of the nuclear luminosities from our photometric decompo- 
sition of the nucleus (§3.3). The above screening process yields 
approximately 200 objects. 

Since the principal aim of this work is to study the intrinsic 
scatter of the MBH~£bui relation, it is imperative that we choose 
objects for which we can obtain the most reliable estimates of 
the two primary quantities of interest, Mbh and L\, u \. Guided 
by this overriding goal, we purposely restrict our sample to the 
upper end of the Mbh distribution. All else being equal, this se- 
lection criterion biases the sample toward more luminous, more 
massive, earlier-type hosts for which we can derive more re- 
liable bulge parameters because the structural decomposition 
will be less complicated than in later-type systems. An added 
benefit of this mass selection is that our sample will consist of 
close analogs of higher-redshift quasars. For concreteness, we 
choose sources with Mbh > 10 78 M Q . This limiois admittedly 
somewhat 



8 Our BH masses assume a smaller geometrical factor than that used in Onken et al. (2004), by a factor of 1.8. Thus, Mbh = 10 7 ' 8 Mq on our scale is equivalent 
to Mbh = 10 8 Mq on the scale of Onken et al. We note that that our main conclusions do not rely on the geometrical factor. This issue is discussed in §5.2. 



Table 1. The sample 



Source Name 


z 


D L 


Line 


FWHM 


Ref. 


log ALKSS 


Ref. 


log *4TSf 


log M BH 


Method 


Ta 


ol/i-Edd 


l0g ^cm 


Ref. 






(Mpc) 




(km s- 1 ) 




(ergs s 1 ) 




(ergs s 


- 1 ) 


(Mq) 








(W Hz" 1 ) 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




(10) 


(11) 




(12) 


(13) 


(14) 



3C 59 
E1821+643 
Fairall 9 

[HB89] 0316-346 
[HB89] 2201+315 
HE 0306-3301 
HE 0354-5500 
HE 1434-1600 
MC 1635+119 
MRK 1048 
MS 0244.8+1928 
MS 0754.6+3928 
MS 1059.0+7302 
MS 1545.3+0305 
OX 169 
PG 0052+251 
PG 0804+761 
PG 0923+201 
PG 0953+414 
PG 1004+130 
PG 1012+008 
PG 1116+215 
PG 1202+281 
PG 1211+143 
PG 1226+023 
PG 1302-102 
PG 1307+085 
PG 1309+355 
PG 1351+640 
PG 1411+442 
PG 1416-129 
PG 1426+015 
PG 1444+407 



0.109 
0.297 
0.047 
0.265 
0.295 
0.247 
0.267 
0.142 
0.146 
0.043 
0.176 
0.096 
0.089 
0.098 
0.211 
0.155 
0.100 
0.190 
0.234 
0.240 
0.186 
0.176 
0.165 
0.080 
0.158 
0.278 
0.155 
0.184 
0.088 
0.089 
0.129 
0.086 
0.267 



499.4 
1526.8 
205.9 
1338.5 
1514.9 
1234.8 
1350.1 
665.3 
685.8 
188.2 
842.9 
435.8 
402.0 
445.5 
1032.6 
732.4 
455.3 
918.0 
1161.0 
1195.0 
896.4 
842.9 
784.7 
359.1 
748.0 
1414.4 
732.4 
885.7 
397.2 
402.0 
599.1 
387.7 
1350.1 



IIq 
H/3 
H/3 

H/3 
H/3 
H/3 
H/3 
H/3 
IT/3 
H/3 
IT/3 
H/3 
H/3 
H/3 
H/3 
H/3 
H/3 

n/3 

II,<3 
IT/3 
H/3 
IT/3 
H/3 
H/3 
H/3 
H/3 
H/3 
H/3 
IIS 
118 
IT/3 
IIS 
TI/3 



9800 
5190 
5900 
7048 
3320 
2346 
2264 
7068 
5271 
7498 
3700 
2123 
7690 
5868 
5030 
4165 
2012 
7233 
3140 
6215 
2604 
3054 
4931 
1317 
2598 
3758 
5058 
3974 
5579 
2398 
6283 
6323 
2944 



43.94 
44.25 



44.06 
43.85 

44.37 



45.18 
44.96 
44.94 
44.88 
45.22 
45.29 
44.47 
45.40 
44.60 
44.75 
45.96 
45.83 
44.88 
45.01 
44.05 
44.63 
45.14 
44.72 
45.20 



44.13 
45.75 
42.98 
44.72 
45.24 
44.79 
44.72 
44.38 
43.50 
43.73 
44.36 
44.71 
43.48 
43.22 
44.72 
44.64 
44.96 
44.74 
45.09 
45.18 
44.44 
45.00 
44.49 
44.58 
45.79 
45.34 
44.64 
44.85 
44.47 
44.50 
44.14 
44.46 
45.00 



8.85 
9.08 
8.15 
8.81 
8.43 
7.89 
7.83 
8.64 
7.92 
8.35 
8.06 
7.76 
8.24 
7.87 
8.52 
8.31 
8.58 
8.84 
8.18 
8.94 
7.80 
8.23 
8.38 
7.90 
8.69 
8.58 
8.38 
8.38 
8.48 
8.39 
8.41 
8.85 
8.21 



+0.20 
-0.20 
+0.20 

0.20 
+0.09 

0.09 
+0.30 

0.30 
+0.30 

0.30 
+0.20 
-0.20 
+0.20 

0.20 
+0.20 
-0.20 
+0.20 

0.20 
+0.20 
-0.20 
+0.20 

0.20 
+0.20 
-0.20 
+0.20 

0.20 
+0.20 

0.20 
+0.20 

0.20 
+0.10 

0.10 
+0.0.1 
-O.O.i 
+0.20 

0.20 
+0.08 
-0.08 
+0.20 

0.20 
+0.20 
-0.20 
+0.20 

0.20 
+0.30 

0.30 
+0.11 

0.11 
+0.08 

0.08 
+0.30 

0.30 
+0.11 

0.11 
+0.30 
-0.30 
+0.20 

0.20 
+0.12 
-0.12 
+0.30 

0.30 
+0.12 
-0.12 
+0.20 

0.20 



S 

s 

R 

S 

s 
s 
s 
s 
s 
s 
s 
s 
s 
s 
s 

K 
R 
S 
R 
S 
S 
S 
S 
R 
R 

s 

R 
S 
S 
R 
S 
R 
S 



0.01 
0.34 
0.09 
0.06 
0.47 
0.57 
0.56 
0.04 
0.05 
0.02 
0.14 
0.73 
0.01 
0.02 
0.19 
0.32 
0.16 
0.07 
0.77 
0.14 
0.32 
0.65 
0.10 
0.50 
1 .32 
0.40 
0.22 
0.25 
0.04 
0.12 
0.12 
0.05 
0.56 



25.27 



26.93 

23.92 
23.30 
24.43 
21.92 

23.45 



26.11 
22.65 
22.75 
22.38 
23.46 
22.85 
22.96 
23.38 
22.79 
24.36 
27.37 
26.27 
22.32 
24.70 
21.99 
22.07 
23.19 
22.31 
22.54 



FT 
X 

2 



a 
a? 



K 

o 
p 



Table 1 — Continued 



Source Name 



z D L 

(Mpc) 
(1) (2) (3) 



Line 



(4) 



FWHM 

(kms- 1 ) 

(5) 



Ref. 



(6) 



(ergs s 1 ) 
(7) 



(8) 



log AL™ a ge logMsH Method L bol /L Edd log P*^ Ref. 
(ergs s- 1 ) (M ) (W Hz" 1 ) 

(9) (10) (11) (12) (13) (14) 



PG 1613+658 


0.129 


599.1 


H/3 


7897 


3 


44.98 


2 


44.80 


8 19+ al ° 

0-1 -0.10 


Ti 


0.44 


23.11 


7 


PG 1617+175 


0.112 


514.2 


H/3 


4718 


3 


44.48 


2 


44.37 


8-siiS:S 


Ti 


0.07 


22.51 


7 


PG 1700+518 


0.292 


1497.0 


IT/3 


1846 


3 


45.63 


2 


45.37 


8.63±°;» 


R 


0.71 


24.29 


7 


PG 2130+099 


0.062 


274.6 


11,9 


2912 


3 


44.46 


2 


44.48 


8 - 4U -0.05 


Ti 


0.08 


22.24 


7 


PHL 909 


0.171 


816.4 


TI/3 


9660 


6 


44.47 


1 


44.78 


9-i2±S:IS 


S 


0.02 


23.02 


8 


PHL 1093 


0.260 


1309.5 


H/3 


7610 


6 


44.30 


1 


44.31 


8-66±g;»g 


S 


0.03 


26.07 


4 


PKS 0736+01 


0.191 


923.4 


H/3 


3259 


4 


44.58 


1 


44.61 


8.08iS;» 


S 


0.23 


26.22 


2 


PKS 1020-103 


0.197 


955.9 


H/3 


7950 


6 


43.88 


1 


44.29 


8-69±g;»g 


S 


0.02 


25.68 


9 


PKS 1217+02 


0.239 


1189.3 


H/3 


3830 


6 


44.89 


1 


44.66 


o nr+0.20 
°- zo -0.20 


S 


0.24 


25.87 


8 


PKS 2135-14 


0.200 


972.3 


H/3 


9294 


4 


45.18 


1 


44.54 


8-96±g;»g 


S 


0.06 


26.16 


9 


PKS 2349-01 


0.174 


832.3 


H/3 


5511 


4 


44.97 


1 


44.63 


8 55+ a2 ° 


s 


0.12 


25.73 


10 


PKS 2355-082 


0.211 


1032.6 


H/3 


7510 


6 


43.98 


1 


44.12 


8 - 56 -0.20 


s 


0.02 


25.30 


10 



Note. — Col. (1): Object name. Col. (2): Redshift. Col. (3): Luminosity distance. Col. (4): Broad emission line used to estimate black hole mass. Col. (5): 
Full width at half maximum of broad emission line. Col. (6): Reference for FWHM: (1) Eracleous & Halpern 1994; (2) Zheng et al. 2002; (3) Peterson et al. 2004; 
(4) Marziani et al. 2003; (5) Letawe et al. 2008; (6) McLure & Dunlop 2001; (7) Boroson 2005. Col. (7): Continuum luminosity at 5100 A derived from observed 
spectrum. Col. (8): Reference for AL^Jj: (1) McLure & Dunlop 2001; (2) Peterson et al. 2004; (3) Grupe et al. 2004; (4) Vestergaard & Peterson 2006. Col. (9): 
Continuum luminosity at 5100 A derived from the best-fit GALFIT model. Col. (10): Black hole mass. Col. (11): Method for estimating Mbh : S — virial method 
based on the single-epoch spectrum; R — reverberation mapping method. Col. (12): Eddington ratio, where Lbol — 9AL>,(5100 A). Col. (13): Radio power at 6 cm; 
for observations done at other wavelengths, we extrapolate to 6 cm assuming /„ <x y~ 05 . Col. (14): Reference for -Pg^: (1) Kellermann & Pauliny-Toth 1973; (2) 
Becker et al. 1991; (3) Mauch et al. 2003; (4) Condon et al. 1998; (5) Barvainis & Antonucci 1989; (6) Becker et al. 1995; (7) Kellermann et al. 1989.; (8) Bennett et 
al. 1986; (9) Griffith et al. 1994; (10) Griffith et al. 1995. 



2 



Table 2. Fitted Parameters 





Detector/Filter 


Exp. 


ffinuc 




Bu 


ge 




Disk 


or Tidal Feature** 




Host 




ai 




Name 


™bul 


Re 


n 


/'f 


"idisk 


Re 


n 


lie 


™host 


W^aper 


B/T 


Qual. 






(s) 


(mag) 


(mag) 


(") 






(mag) 


(") 






(mag) 


(mag) 








(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


3C 59 


WF2/F675W 


1800 


15.94 


15.87 


3.41 


4.00 


21.93 


< 20.5 








15.87 


15.9 


1(> 0.99) 


0.01 


X 


E1821+643 


WF3/F791W 


3900 


13.51 


15.97 


4.31 


4.00 


22.53 


< 19.0 










15.97 


15.6-16.3 


1(> 0.94) 


0.17 


X 


Fairall 9 


HRC/F550M 


1020 


17.97 


15.83 


0.04 


4.00 


12.18 


16.45 


1.78 


1.00 


20 


.40 


15.35 


15.2-15.5 


0.64 


0.03 


X 


[HB89] 0316-346 


WF3/F606W 


1400 


16.53 


18.57 


1.34 


3.00 


22.44 


19.92 


1.28 


0.73 


23 


.01 


18.06 




0.63 


0.18 


s 


[HB89] 2201+315 


PC/F702W 


560 


15.41 


17.34 


0.78 


4.00 


20.18 


< 20.0 










17.34 




1(> 0.92) 


0.09 


s 


[HB89] 2201+315 


PC/F555W 


2700 


15.95 


17.87 


1.38 


4.00 


21.96 


< 20.0 










17.87 




1(> 0.88) 


0.01 


s 


HE 0306-3301 


HRC/F606W 


990 


16.29 


19.49 


0.68 


4.00 


22.04 


18.47 


1.33 


1.00 


21 


.79 


18.11 


17.5-18.2 


0.28 


0.06 


X 


HE 0354-5500 


HRC/F606W 


1200 


16.56 


18.19 


0.58 


4.00 


22.40 


20.98 










18.11 


17.4-18.2 


0.93 


0.33 


X 


HE 1434-1600 


HRC/F606W 


990 


16.43 


17.54 


1 .85 


4.00 


22.27 


< 21.0 










17.54 


17.2-17.4 


1(> 0.96) 


0.01 


X 


MC 1635+119 


WF2/F675W 


1800 


18.19 


16.80 


1 .76 


4.00 


21.42 


< 21.0 










16.80 


16.7 


1(> 0.98) 


0.01 


X 


MRK 1048 


PC6/F785LP 


230 


14.91 


13.69 


1 .61 


4.00 


18.12 


14.52 


14.05 


1.00 


22 


.95 


13.28 


13.3 


0.68 


1 .86 


X 


MS 0244.8+1928 


WF2/F675W 


1800 


16.82 


17.55 


2.07 


4.00 


22.52 


< 20.5 










17.55 


17.5 


1(> 0.94) 


0.04 


X 


MS 0754.6+3928 


PC/F814W 


600 


13.70 


15.94 


2.34 


4.00 


21.18 


< 18.0 










15.94 


15.9 


1(> 0.87) 


0.09 


X 


MS 1059.0+7302 


PC/F814W 


600 


16.61 


17.02 


0.38 


3.94 


18.33 


15.89 


2.76 


1.00 


20 


.80 


15.57 


15.6 


0.26 


0.00 


X 


MS 1545.3+0305 


PC/F814W 


600 


17.57 


16.30 


1.09 


4.00 


19.89 


16.20 


3.23 


1.00 


21 


.44 


15.50 


15.5 


0.47 


0.36 


X 


OX 169 


WF2/F675W 


1800 


16.02 


17.33 


1 .05 


4.00 


20.82 


19.05 










17.13 


16.8-16.9 


0.83 


0.29 


X 


PG 0052+251 


WF2/F675W 


1800 


15.48 


16.95 


1 .38 


4.00 


21.02 


17.72 


3.99 


0.20 


22 


.86 


16.51 


16.3-16.4 


0.67 


0.22 


X 


PG 0804+761 


WFC/F625W 


1590 


13.95 


16.83 


1.75 


4.00 


21.44 


19.06 


2.17 


0.22 


22 


.88 


16.70 


16.4 


0.89 


0.01 


X 


PG 0923+201 


WF2/F675W 


1800 


15.66 


17.22 


2.04 


4.00 


22.16 


< 20.0 










17.22 


16.8-17.2 


1(> 0.93) 


0.00 


X 


PG 0953+414 


WF2/F675W 


1800 


15.11 


18.03 


2.61 


4.00 


23.50 


< 21.5 










18.03 


17.9-18.7 


1(> 0.96) 


0.12 


X 


PG 1004+130 


WF2/F675W 


1800 


15.04 


16.79 


2.05 


4.00 


21.73 


< 21.5 










16.79 


16.7-16.9 


1(> 0.99) 


0.04 


X 


PG 1012+008 


WF2/F675W 


1800 


16.33 


17.08 


2.96 


3.24 


22.72 












17.08 


15.9-16.0 


1 


0.31 


X 


PG 1116+215 


WF3/F606W 


2200 


15.27 


17.01 


1.99 


4.00 


21.90 


< 20.5 










17.01 


16.6-16.9 


1(> 0.96) 


0.02 


X 


PG 1202+281 


WF3/F606W 


1800 


16.25 


17.63 


1 .33 


4.00 


21.64 


< 20.5 










17.63 




1(> 0.93) 


0.00 


X 


PG 1211+143 


WFC/F625W 


1590 


14.33 


16.94 


0.95 


5.84 


20.42 


18.21 
19.60 


4.72 
1.40 


0.32 
0.15 


23 
23 


.81 
.43 


16.57 


15.7-16.0 


0.72 


0.03 


X 


PG 1226+023 


WF3/F606W 


2200 


12.86 


16.06 


6.28 


4.00 


23.44 












16.06 




1 


0.13 


p 


PG 1302-102 


WF3/F606W 


1800 


15.12 


17.73 


1 .03 


4.00 


21.19 


< 21.0 










17.73 




1(> 0.95) 


0.00 


s 


PG 1302-102 


PC/F702W 


1800 


14.94 


17.34 


2.03 


4.00 


22.27 


< 20.5 










17.34 




1(> 0.95) 


0.01 


s 


PG 1307+085 


WF3/F606W 


1400 


15.72 


17.72 


0.94 


4.00 


20.96 


< 21.5 










17.72 


17.7 


1(> 0.97) 


0.07 


s 


PG 1309+355 


WF3/F606W 


2100 


15.50 


16.74 


1.58 


4.00 


21.12 


< 21.0 










16.74 




1(> 0.98) 


0.01 


s 



FT 
X 

en 

2 



a 
a? 



o 




Table 2 — Continued 





Detector/Filter 


Exp. 


?Tlnuc 




Bu 


ge 




Disk 


or Tidal Feature a 




Host 




ai 




Name 


"Ihul 


Re 


n 


lie 


™disk 


Re 


n 


/'f 


™host 


diaper 


BjT 


Qual. 






(<0 


(mag) 


(mag) 


(") 






(mag) 


(") 






(mag) 


(mag) 








(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


PG 1351+640 


WFC/F625W 


1800 


14.85 


16.07 


1.96 


3.81 


20.90 


18.75 


1.41 


1.00 


22.19 


15.98 


15.5 


0.92 


0.09 


N 


PG 1411+442 


WFC/F625W 


1680 


14.76 


16.96 


1.12 


4.00 


20.60 


16.68 


8.08 


1 .38 


24.06 


16.06 


15.8-16.4 


0.43 


0.21 


N 


PG 1416-129 


PC/F814W 


600 


15.84 


17.01 


3.09 


4.00 


22.84 


< 20.0 








17.01 


16.1-16.8 


1(> 0.94) 


0.00 


S 


PG 1426+015 


PC/F814W 


600 


14.02 


16.08 


0.61 


2.12 


18.08 


15.26 


4.11 


0.71 


20.87 


14.84 


14.3-14.7 


0.32 


0.12 


N 


PG 1444+407 


WF3/F606W 


2300 


15.81 


18.68 


0.85 


4.00 


21.71 


19.08 


2.47 


0.31 


23.26 


18.10 


17.5-17.8 


0.59 


0.10 


N 


PG 1613+658 


WFC/F625W 


1860 


14.87 


15.92 


3.00 


4.00 


21.69 


16.25 


4.89 


1.00 


22.39 


15.31 


15.0-15.1 


0.58 


0.10 


N 


PG 1617+175 


WFC/F625W 


1590 


15.69 


17.70 


1 .37 


4.00 


21.78 


19.51 


1 .83 


0.22 


22.98 


17.51 


17.4-17.5 


0.84 


0.10 


N 


PG 1700+518 


PC/F547M 


2000 


15.45 


17.78 


1.15 


4.00 


21.47 


19.82 








17.63 


16.0-17.3 


0.87 


0.39 


N 


PG 2130+099 


WF2/F555W 


1800 


14.47 


17.34 


2.31 


0.44 


21.51 


16.49 


12.85 


0.33 


24.28 


16.08 


16.0 


0.31 


0.23 


S 


PHL 909 


WF2/F675W 


1800 


15.51 


16.60 


2.72 


4.00 


22.16 


< 21.5 








16.60 


16.1-16.6 


1(> 0.99) 


0.10 


N 


PHL 1093 


WF2/F675W 


1800 


17.27 


17.23 


2.24 


4.00 


22.37 


< 19.5 








17.23 


17.0-17.2 


1(> 0.89) 


0.07 


N 


PKS 0736+01 


WF2/F675W 


1800 


16.20 


16.93 


2.77 


4.00 


22.54 


< 20.5 








16.93 


16.7-17.0 


1(> 0.96) 


0.04 


N 


PKS 1020-103 


WF2/F675W 


1800 


16.82 


17.18 


1 .68 


4.00 


21.69 


< 22.0 








17.18 


17.1-17.2 


1(> 0.99) 


0.05 


N 


PKS 1217+02 


WF2/F675W 


1800 


16.31 


17.28 


2.21 


4.00 


22.39 


< 20.5 








17.28 


17.3 


1(> 0.95) 


0.03 


X 


PKS 2135-14 


WF2/F675W 


1800 


16.27 


17.14 


2.97 


4.00 


22.90 


< 21.0 








17.14 


16.9-17.0 


1(> 0.97) 


0.04 


X 


PKS 2349-01 


WF2/F675W 


1800 


16.02 


16.41 


3.04 


4.00 


22.22 


17.12 


3.30 


0.88 


22.34 


15.96 


16.0 


0.65 


0.69 


X 


PKS 2355-082 


WF2/F675W 


1800 


17.33 


17.05 


2.13 


4.00 


22.09 


< 20.5 








17.05 


16.9-17.1 


1(> 0.96) 


0.05 


X 



2 



a We list R e , n, and ji e only for unambiguous disks; otherwise, we attribute the component to a tidal features, for which only their apparent magnitude is given. 



Note. — Col. (1): Object name. Col. (2): HST detector and filter. Col. (3): Exposure time. Col. (4): Apparent nuclear magnitude at the observed filter. 
Col. (5): Apparent bulge magnitude at the observed filter. Col. (6): Effective radius of bulge. Col. (7): Sersic index for bulge. Col. (8): Surface brightness of 
bulge at the effective radius (mag arcsec -2 ). Col. (9): Apparent magnitude of disk or tidal feature. Col. (10): Effective radius of disk. Col. (11): Sersic index for 
disk. Col. (12): Surface brightness of disk at the effective radius (mag arcsec -2 ). Col. (13): Apparent host magnitude from combining the luminosity of the bulge 
(Col. 5) and of the disk or tidal features (Col. 9). Col. (14): Apparent host magnitude from aperture photometry after subtracting the PSF, determined either by 
the best-fit PSF or by scaling up to the maximum value of image. Col. (15): Bulge-to-total luminosity ratio of host galaxy; upper limit, if measurable, is shown 
in parentheses for elliptical galaxies (see text for details). Col. (16): Amplitude of the first Fourier mode. Col. (17): Quality of image: N = normal image; S = 
saturated image; P— poor image (pixels near object are useless due to bleeding columns). 



Black Hole Mass vs. Bulge Luminosity Relation 



Table 3. Derived Properties 



Name 


Av 


-Mfl.nuc 


M H ,bul 


-Mb, host 


(i) 


(mag) 
(2) 


(mag) 
(3) 


(mag) 
(4) 


(mag) 
(5) 


3C 59 

E1821+643 


0.211 

0.141 


-22.71 
-27.15 


-22.66l°;3 
— 24.63+g'g 


— 22.66±55:l 

-24.63+55-3 


Fairall 9 

[HB89] 0316-346 
[HB89] 2201+315 
[HB89] 2201+315 
HE 0306-3301 


0.089 
0.056 
0.410 
0.410 
0.046 


-18.98 
-24.50 
-25.87 
-25.91 
-24.56 


-21.16+™ 
— 22.86±°;| 
-23.80+™ 
-24.571™ 
-21.Tl±g"i 


-21.561™ 
-23.141™ 
-23.801™ 
-24.571™ 
-22.651™ 


HE 0354-5500 


0.053 


-24.49 


-23.261™ 


— 23.30±°;1 


HE 1434-1600 
MC 1635+119 


0.417 
0.171 


-23.38 
-21.27 


-22.09l°;3 

-22.47 + 55-3 


— 22.o9±55:l 

-22.47155-3 


MRK 1048 

MS 0244.8+1928 


0.110 
0.368 


-21.37 
-23.26 


-21.94+™ 
-22.21±°;1 


— 22.49+55 1 

-22.21!™ 


MS 0754.6+3928 


0.220 


-24.22 


-21.92±°;1 


-21.92!™ 


MS 1059.0+7302 


0.173 


-21.12 


-20.65+™ 


-22.10!™ 


MS 1545.3+0305 


0.377 


-20.48 


-21.61+™ 


-22.40!™ 


OX 169 

PG 0052+251 

PG 0804+761 


0.367 
0.157 
0.115 


-24.50 
-24.12 
-24.68 


-22.93+™ 
-22.471™ 
-21.66+°-* 


— 23.io±55:l 
-22.89155-3 

-21.781™ 


PG 0923+201 


0.141 


-24.44 


-22.76J:°-3 


-22.76155-3 


PG 0953+414 


0.042 


-25.41 


-22.S3±°;| 


-22.531™ 


PG 1004+130 


0.127 


-25.61 


-23.841™ 


-23.841™ 


PG 1012+008 


0.116 


-23.71 


-22.83l°;3 


-22.83l5J;3 


PG 1116+215 


0.075 


-24.78 


-23.20l°;3 


— 23.2o±55:l 


PG 1202+281 


0.070 


-23.64 


-22.40+55-5 


-22.40155-5 


PG 1211 + 143 


0.115 


-23.76 


-21.001551 


-21.361™ 


PG 1226+023 


0.068 


-26.92 


-23.85±!5;1 


-23.851™ 


PG 1302-102 
PG 1302-102 


0.141 
0.141 


-26.10 
-26.00 


-23.85+°; = 

-23.63+5J-5 


-23.851™ 
-23.63+5J-5 


PG 1307+085 


0.112 


-24.05 


-22.14+™ 


-22.141™ 


PG 1309+355 


0.040 


-24.63 


-23.591™ 


-23.591™ 


PG 1351+640 


0.068 


-23.43 


-22.1llg1 


-22.191™ 


PG 1411+442 


0.028 


-23.52 


-21.241551 


-22.111™ 


PG 1416-129 
PG 1426+015 


0.311 
0.105 


-22.82 
-23.61 


-21.571™ 
-21.5llg1 


-21.571™ 
-22.751™ 


PG 1444+407 


0.046 


-25.23 


-22.77+55-4 


-23.101™ 


PG 1613+658 


0.088 


-24.33 


-23.211551 


-23.73155-3 


PG 1617+175 


0.138 


-23.22 


-21.081551 


-21.241™ 


PG 1700+518 


0.116 


-25.97 


-24.381551 


-24.441™ 


PG 2130+099 


0.147 


-23.44 


— 20.68±55;1 


-21.611™ 


PHL 909 
PHL 1093 
PKS 0736+01 


0.151 
0.096 
0.422 


-24.34 
-23.53 
-24.12 


— 23.08 + 55:l 

— 23.62±° 1 

-23.051™ 


— 23.08ioJ 
— 23.62±° | 

"23.051™ 


PKS 1020-103 


0.153 


-23.38 


-22.8915J-3 


— 22.89±°'.l 


PKS 1217+02 


0.070 


-24.29 


-23.341™ 


— 23.34±°:1 



KIM et al. 



Table 3 — Continued 



Name 

(i) 




Av 

(mag) 

(2) 


(mag) 
(3) 


A^B.bul 

(mag) 
(4) 


M R,hoat 

(mag) 
(5) 


PKS 2135- 
PKS 2349- 
PKS 2355- 


-14 
01 
-082 


0.169 
0.091 
0.133 


-23.98 
-23.81 
-23.02 


-23.05 + o| 
-23.32l°i 
-23.21±g;| 


-23.051°;^ 
-23.73±o | 
-23.2lj;°j 



Note. — Col. (1): Object name. Col. (2): Galactic extinction in 
the V band. Col. (3): .R-band absolute magnitude of the nucleus 
after correcting for Galactic extinction. Col. (4): .R-band absolute 
magnitude of the bulge and its error (see text for details). Col. (5): 
.R-band absolute magnitude of the entire host. 



Black Hole Mass vs. Bulge Luminosity Relation 



Table 4. log M BH /M = a - 0.50M R , 



bill 



Sample 


a 


Aa a 


Scatter 


All 

High L bol /L E dd (> 0.1) 
Low L bol /L Edd (< 0.1) 


-2.87 ±0.04 

-3.08 ± 0.05 
-2.53 ±0.06 


-0.27 

-0.48 
0.07 


0.41 

0.31 
0.30 



Elliptical (B/T = 1) 
Bulge-dominated (B/T > 0.5) 
Disk-dominated (B/T < 0.5) 



-3.04 ±0.05 -0.44 0.27 
-2.74 ±0.07 -0.14 0.43 
-2.26 ±0.11 0.34 0.37 



No Interaction (ai < 0.3) -2.78 ± 0.04 -0.18 0.37 
Interaction (ai > 0.3) -3.33 ±0.10 -0.73 0.32 



Radio-loud (R > 10) 
Radio-quiet (R < 10) 



-3.05 ±0.07 -0.45 0.22 
-2.78 ±0.05 -0.18 0.45 



z< 0.15 
z> 0.15 



-2.42 ±0.04 0.18 0.44 
-3.10 ±0.04 -0.50 0.28 



a Ao; compared to the relation of normal galaxies (a = —2.60). 
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c 

a 



X 




FIG. 3. — Distribution of the differences in host galaxy magnitudes from 
the 2-D imaging-fitting simulations of Kim et al. (2008). Artificial images of 
AGN host galaxies with a range of input parameters were generated, and GAL- 
FIT was used to recover the input parameters. The open histograms denote the 
eiTors for idealized conditions in which the fitting was done with the same PSF 
as used for making the input images. The hatched histograms give the errors 
for the more realistic situation that accounts for PSF mismatch (see Kim et al. 
2008 for details). PSF mismatch causes significant systematic errors. 

arbitrary, but it yields a sizable sample of 45 objects. Our ob- 
jects (Table 1) span ~ 2 orders of magnitude in accretion rate 
(Eddington ratio) over a relatively narrow range (^ 1 dex) in 
Mbh (Fig. 1) and redshift (Fig. 2). From radio data assembled 
from the literature, 33% (15/45) of the sample is radio-loud, de- 
fined by R > 10, where R = fJ6 cm)// 1/ (4400 A) (Kellermann 
et al. 1989). 

3. ANALYSIS 
3.1. Structural Decomposition 



outweighs most other concerns, including the choice of actu- 
ally observed (stellar) or synthetic (Tiny Tim; Krist 1995) PSFs. 
They demonstrate that the undersampling problem can be sig- 
nificantly alleviated by broadening both the science image and 
the PSF image to critical sampling [full width at half maximum 
(FWHM) ~ 2 pixels]. This is the strategy we adopt here. Al- 
though the choice of real versus synthetic PSFs is secondary, 
we use observed PSF stars whenever possible. When these are 
not available, we use Tiny Tim PSFs. 

As in our companion paper, we use an updated version of 
GALFIT (C. Y. Peng et al., in preparation^ The code simul- 
taneously fits multiple components to model the host galaxy, 
with the freedom to use Fourier modes to accommodate com- 
plex, nonaxisymmetric features such as tidal distortions or even 
spiral arms. These improvements allow us to obtain a more ac- 
curate decomposition of the structural components of the host, 
an important consideration for our aim of deriving robust bulge 
luminosities. 

We model the active nucleus with a synthetic PSF and the 
host galaxy with ellipses represented by a Sersic (1968) func- 
tion: 



I{r) = I e exp 



-b„ 



l/n 



(1) 



where r e is the effective radius, I e is the intensity at r e , n is the 
Sersic index, and b„ satisfies 



I(r)2irrdr = 2 / I(r)2Trrdr. 



Our companion paper (Kim et al. 2008) discusses in detail 
our technique for decomposing the HST images of AGN host 
galaxies. We performed extensive simulations to quantify the 
performance of the two-dimensional (2-D) image-fitting code 
GALFIT (Peng et al. 2002) under conditions typically encoun- 
tered in AGN host galaxy images contained in the HST archive, 
similar to those analyzed in this study. We paid particular atten- 
tion to quantifying systematic uncertainties in estimating the 
photometric parameters of the bulge component and devised 
strategies for mitigating them. 

In the regime where the bright, unresolved active nucleus 
dominates over the light of the host galaxy, our ability to ex- 
tract the bulge luminosity depends sensitively on the properties 
of the point-spread function (PSF). PSF mismatch systemati- 
cally biases the derived bulge luminosities to high values, by as 
much as 0.5-1 mag (Fig. 3). PSF mismatch occurs as a result of 
variations in time, location on the detector, and, to a lesser ex- 
tent, differences in the spectral energy distribution between the 
PSF star and the science target. The dominant effect, however, 
comes from the fact that HST PSFs are undersampled. Because 
the PSFs are not Nyquist-sampled, it is impossible to shift by 
a sub-pixel unit to perfectly align the PSF star with the AGN 
core. Kim et al. (2008) show that this fundamental problem 

^http : //www . ociw . edu/~ peng/work/ga lf it /newf eatures ■ html| 



(2) 



The Sersic function reduces to an exponential profile for n = 1 
and a de Vaucouleurs (1948) profile for n = 4. 

Several different fits are done for each object. We determine 
the sky value from a growth curve analysis of the image. We 
begin with the simplest possible option of modeling the host 
with a single Sersic component, allowing the index n to be free 
as well as fixing it to n = 1 and n = 4. Our simulations show 
that when the AGN is much brighter than the bulge allowing 
n to be free can sometimes lead to spurious results. In these 
situations, it is better to fix n to specific values in the fit and 
then empirically bracket the resulting host luminosity from the 
range of acceptable models, as judged from the relative change 
of x 2 and visual examination of the model residuals. With the 
Sersic index fixed, the program solves for the following free 
parameters: for the host galaxy, these are the position, effective 
radius R e , surface brightness /i e at R e , axis ratio, and position 
angle; for the AGN, these are the position and nuclear magni- 
tude, m nuc . Figure 56 gives an example of a single-component 
fit for PKS 2355-082. 

We perform two-component (bulge+disk) fits to objects that 
clearly have a disk component in the original image or that 
show significant extended structure in the residual image after 
subtraction of a single-component model. We try three cases: 
(1) n = 4 for bulge and n = 1 for disk; (2) n = free for bulge 
and « = 1 for disk; and (3) n = free for bulge and n = free for 
disk. We permit the bulge to have n < 4 to allow for the possi- 
bility of pseudobulges (Kormendy & Kennicutt 2004), and the 
disk can deviate from a pure exponential profile if it is severely 
distorted, if it has bar or ring-like structures, or if it represents 
tidal features and arcs. In most cases when a disk is present, it 
tends to have a smaller Sersic n, a larger R e , and often a smaller 
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FIG. 4. — Examples of non-interacting and interacting objects, based on the strength of the first Fourier mode (ay) measured with GALFIT. From left to right, 
the six objects are arranged in the order of increasing ay. The first three, with ay < 0.1, show essentially no signs of perturbation. The last three, with ay > 0.1, 
are increasingly disturbed morphologically. In each column, we show, from top to bottom, the original image, model, and residual. The units of the images are in 
arcseconds. 



axial ratio than the bulge. The disk component adds four addi- 
tional free parameters, namely position, R e , [i e , axis ratio, and 
position angle. Figure 24 gives an example of a bulge+disk fit 
for MS 1059.0+7302. 

The choice between a single-component and a double- 
component fit for the host galaxy is not always clear. Adding 
extra free parameters obviously yields a better fit. The question 
is whether the extra component is clearly required and physi- 
cally meaningful. For example, when PSF mismatch is partic- 
ularly severe, the extra component might simply be attempting 
to account for the large residuals from the poor PSF model. In 
such situations the "extra" component tends to have unusual 
properties such as suspiciously tiny R e or extreme values of n. 
In practice, complicated fits often require some degree of judg- 
ment call, but we have tried to err on the side of caution and 
generally invoke extra components only when they are abso- 
lutely needed. In most cases, we admit an additional compo- 
nent only if it is clearly visible in the original image or in the 
residual image. 

A disk component may be present but undetectable in shal- 
low images (e.g., Bennert et al. 2008). It is thus very useful to 
place an upper limit on the disk component even if no disk is re- 
quired by the best-fit model of the HSr image. For objects with 
no directly detected disks, we derive upper limits for the disk 
component by placing artificial (face-on) disks covering a wide 
range of luminosity on the science image, assuming that the 
host galaxy follows the relation between bulge-to-total light ra- 
tio (B/T) and the ratio of disk size to bulge effective radius de- 
rived from nearby early-type galaxies (de Jong et al. 2004). We 
then fit the simulated images with two-component (bulge+disk) 



models. The disk luminosity at which the program fails to re- 
cover the input value gives an estimate of the upper limit for the 
disk component. We did not attempt to derive disk upper lim- 
its for systems that are exceptionally complicated (e.g., highly 
distorted, close companions, etc.). 

The updated version of GALFIT also has the ability to model 
spiral arms in the disk. The spiral structure is created by a hy- 
perbolic tangent rotation function with the following parame- 
ters: bar length, outer spiral radius, rotation rate, sky inclina- 
tion, and position angle. The details of the spiral structures are 
created by high-order Fourier modes. Figure 41 shows the fit 
for PG 141 1+442, whose disk component shows two prominent 
spiral arms. 

A unique aspect of our analysis is that we attempt to quantita- 
tively estimate the degree to which the host galaxy exhibits non- 
axisymmetric distortions. Morphological disturbances may be 
signatures of recent mergers or tidal interactions, which might 
trigger or enhance AGN fueling. While a variety of techniques 
have been devised to characterize morphological asymmetry in 
inactive galaxies (e.g., Conselice et al. 2000; Lotz et al. 2004), 
they cannot be readily extended to galaxies containing bright 
AGNs because the central point source can dominate the asym- 
metry signal. The latest version of GALFIT implements asym- 
metry parameters as an integral part of the image-fitting pro- 
cess. This is accomplished by introducing higher-order Fourier 
modes to change the shape of the host galaxy from axisym- 
metric ellipses into more complicated shapes. All the while 
the light profile of the host galaxy model would still decline 
as a Sersic profile in every direction from the peak. In this 
scheme, the strength of an external perturbation on the host 
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galaxy would sensitively register as high-amplitude Fourier 
modes, with phase angles that reflect the direction of the per- 
turbation. 

If the residual image shows significant nonaxisymmetry, we 
adopted a Fourier component to fit it. The Fourier mode has the 
following form: 



r(x,y) = r (x,y) 



1 + / a m cos(m(6 + 7r m )) 



in- 1 



(3) 



In this expression, 9 = arctan[(y-y c )/(x-x c )g], where (x c ,y c ) 
is the centroid of the ellipse, q is the axis ratio, ro(x,y) is the 
generalized ellipse, a m is the amplitude for mode m, and ir m is 
the phase angle for mode m. The Fourier mode is always cou- 
pled with a general single (e.g., Sersic) component and shows 
how much a component is perturbed from the perfect ellip- 
soid. Thus, the Fourier mode allows us to quantify the de- 
gree of asymmetry. In principle, we can use an infinite num- 
ber of Fourier modes, but in practice we find that four modes 
(m = 1,3,4,5) are enough to fit the asymmetrical structures en- 
countered in our sample. Figure 4 illustrates a series of objects 
with increasing strengths in a\. Sources with a\ <; 0.1 show lit- 
tle to no obvious signs of morphological perturbation, whereas 
those with a\ ;> 0.1 appear increasingly disturbed. 

Despite the significant new features of the updated version 
of GALFIT, we note that the derived values of many standard 
parameters (e.g., size and luminosity) are not substantially dif- 
ferent between the original and new versions of the code. In 
situations where there are differences, the new features of the 
code allow for better convergence, especially when it comes 
to multi-component decompositions. The updated version of 
GALFIT has been tested extensively by us, but for the sake of 
brevity we defer a full discussion of the technical details to an 
upcoming paper (C. Y. Peng et al., in preparation). 

Table 2 summarizes the results of the structural decomposi- 
tion. For each object, we list the best-fit nuclear and photomet- 
ric parameters for the bulge, disk or tidal feature, and for the 
overall host galaxy; the parameter a\ is also tabulated. 

3.2. Bulge Luminosities 

The error bars on the bulge measurements are influenced by 
a number of systematic uncertainties. Using the simulations in 
Kim et al. (2008) as a guide, the final error budget on the bulge 
luminosity was estimated as follows. For sources with L\, u i/L nuc 
> 0.2, ctrs ±0.3 mag, whereas er« ±0.4magifLbui/^nuc <0.2. 
On top of these values, additional uncertainties are introduced 
if bulge-to-disk decomposition is required (^0.1 mag), if satu- 
ration occurs (~ 0.2 mag), or if the image contains substantial 
inner fine structure (~ 0.3 mag). We also determine the error 
of the host luminosity (mh ost ) from the simulations at a given 
Lhost/^nuc- According to Kim et al., the error on the nucleus 
magnitude is ~ ±0.1 mag. 

Because of the complexity of the GALFIT decomposition, it 
is worthwhile to cross-check our 2-D parametric fits with a non- 
parametric estimate of the host luminosity (see, e.g., Greene et 
al. 2008). To perform this test, we remove the nucleus from 
each source simply by subtracting a shifted, scaled PSF model 
from the peak of the AGN core. After masking out obvious 
companions and foreground stars, we sum up the remaining 
flux to estimate the total host galaxy magnitude. As a separate 
test, we compute the host galaxy flux after subtracting the PSF 
component derived from the best-fitting GALFIT model for the 



entire image. These two tests give the range of values tabulated 
as w aper in Table 2. Comparison of these estimates with the host 
magnitudes obtained from the parametric fits (mhost) shows rea- 
sonably good agreement for the majority of the sources. The 
few cases in which m aper is substantially brighter than mh os t can 
be attributed to large residuals from PSF mismatch and contam- 
ination from neighboring sources. 

A significant number of the sources in our sample overlap 
with those studied by Dunlop et al. (2003), affording an inde- 
pendent, external check of our analysis. Dunlop et al. also per- 
formed 2-D decomposition of their sample, but they fitted the 
host galaxies with only a single component, modeled as either a 
classical de Vaucouleurs (n = 4) bulge or an exponential (n = 1) 
disk. After accounting for differences in the adopted cosmo- 
logical parameters, we find, not surprisingly, that for bulge- 
dominated sources our bulge magnitudes generally agree well 
(to within 0. 1 mag) with those given by Dunlop et al. The ex- 
ceptions are objects with large, nearby neighbors and sources 
with multiple components. Whereas we perform a simultane- 
ous 2-D fit of all nearby sources that could potentially affect our 
target of interest, Dunlop et al. simply masked them out. This 
could lead to systematic errors in the derived properties of the 
primary host. 

A particularly striking example is PG 1012+008, which is an 
obviously interacting system consisting of three galaxies. Si- 
multaneously accounting for the subcomponents, including an 
off-centered disk, our best-fit model yields a bulge with mt, u i = 
17.1 (F675W) and an effective radius of R e = 2."96, or 9.3 kpc. 
By contrast, Dunlop et al. obtain m\, u i = 16.4 mag and R e = 
5. "75, which corresponds to 18 kpc using our assumed distance 
of 896 Mpc. For sources that clearly contain both a bulge and 
a disk (e.g., PKS 2349-01), our two-component fits yield more 
robust bulge luminosities. Lastly, for completeness, we note 
that our final nuclear magnitudes (Table 3), converted to the R 
band, are systematically brighter by 0.4 mag compared to those 
given in Dunlop et al. This difference can be traced to the differ- 
ent assumptions used for calculating the ^-correction. Dunlop 
et al. assumed that the AGN spectrum can be represented by 
a single power law (/„ oc v~ 2 ), whereas we use the empirical 
quasar composite spectrum of Vanden Berk et al. (2001). 

Five of our objects (MS 0754.6+3928, MS 1059.0+7302, MS 
1545.3+0305, PG 1416-129, and PG 1426+015) overlap with 
the sample studied by Schade et al. (2000), who also performed 
2-D fits to derive photometric parameters for the host galaxies. 
The two studies show significant difference in the sense that 
Schade et al. tend to underestimate the nuclear magnitudes on 
average by 0.2 mag and to overestimate the bulge magnitudes 
by ~0.7 mag. Some objects show particularly striking disagree- 
ment. In our analysis, the host galaxy of MS 1059.0+7302 is 
well described by a bulge+disk model. Our best fit yields the 
7«bui = 17.02 mag and R e = 0."38 for the bulge and m^^ = 15 .89 
mag and R e = 2."76 for the disk. By contrast, Schade et al. find 
7«bui = 15.53 mag and R e = 3."47 for the bulge and m^^ = 16.68 
mag and R e = 1 "28 for the disk. We attribute the discrepancy 
between our results and those of Schade et al. to a difference 
in methodology. Schade et al. fitted the HST images simultane- 
ously with ground-based images. Although the ground-based 
images are deeper, they have a much broader and less stable 
PSF than the HST images; it is difficult to know how this ef- 
fect impacts the fitting results. Other differences stem from the 
model adopted in the fit. In our work, PG 1426+015 is best fit 
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FIG. 5. — Correlation between BH mass and absolute R-band magnitude 
for the bulge. The BH masses are derived from reverberation mapping (red 
circles) or from single-epoch spectra (blue stars), and the bulge luminosities 
are based on the GALFIT decomposition. The values of Mbh estimated from 
single-epoch spectra are assumed to have a systematic uncertainty of 0.2 dex. 
The best fit is plotted as a solid line, and the intrinsic scatter is denoted by 
the dashed lines. The Mbh — Ami relation for inactive galaxies is shown by the 
thick hatched line. 

with a two-component bulge+disk model with a significant ci\ 
Fourier mode, whereas Schade et al. employed only a single- 
component bulge for the host. If we adopt a single-component 
model, our results agree well with those of Schade et al. 

The bulge magnitudes listed in Table 2 were derived from im- 
ages taken in different filters. For our subsequent analysis, we 
need to convert the magnitudes to a single standard bandpass at 
z = 0. For ease of comparison with the Mbh _ ^bui relation of 
inactive galaxies (Bettoni et al. 2003), we choose the R band as 
the reference. We perform the color conversion of the observed 
magnitude in the various HST filters to the R band and apply 
^-correction using galaxy template spectra from Calzetti et al. 
(1994) and Kinney et al. (1996). We assume that the bulge 
component has the spectrum of an elliptical galaxy and that the 
disk component is approximated by a late-type (Sc) galaxy. For 
the images taken in the F814W filter, we employ the template 
spectrum of a starburst galaxy for the disk component because 
the template spectrum of a late-type galaxy is unavailable in 
this wavelength regime. The final /?-band absolute magnitudes 
are given in Table 3. 

3.3. Black Hole Masses 

The BH masses for type 1 AGNs can be estimated from the 
virial product Mbh ~ fRv 2 /G, where R is the radius of the 
broad-line region (BLR), v is the line width of the BLR gas 
represented by FWHM H /3, the FWHM of the broad H/3 line, 
and / is a factor of order unity that depends on the structure, 
dynamics, and inclination angle of the BLR. Direct measure- 
ments of R through reverberation mapping are available only 
for a small number of sources (Peterson et al. 2004), but for- 
tunately this quantity can be estimated from the correlation be- 
tween R and luminosity (Kaspi et al. 2000, 2005). The virial 
product is, however, uncertain by the normalizing factor /. As- 
suming that the BLR is spherical and has an isotropic velocity 



field, / = 0.75, and the latest radius-luminosity relation from 
Bentz et al. (2006) yields, 



M BH = 5.5 x 10 6 M G 



/ XL 5 



V10 44 ergss- 



0.52 



/FWHM 



Hfi 



V10 3 kms- 



2.0 



(4) 

where AL5100 is the continuum luminosity at 5100 A (see 
Greene & Ho 2007b for details). 

The continuum luminosity can be estimated either through 
spectrophotometry or through our image analysis. Spectropho- 
tometry has the advantage that the specific continuum flux at 
the desired wavelength can be directly measured without mak- 
ing assumptions about the spectral shape. On the other hand, 
accurate absolute spectrophotometry is nontrivial to achieve 
and is rarely available for most objects in the literature. More- 
over, ground-based apertures invariably blend the nucleus with 
at least part of the host. By contrast, our careful image de- 
composition yields a clean, unambiguous measurement of the 
nuclear continuum. Our nuclear magnitudes have a typical un- 
certainty (dominated by systematic effects from PSF mismatch) 
of ~0. 1 mag. We need to assume a spectrum (we choose the 
quasar template from Vanden Berk et al. 2001) in order to es- 
timate the continuum luminosity at 5100 A, but the amount of 
extrapolation for our filters is minimal. A more significant un- 
certainty comes from temporal variability between our photo- 
metric measurements and the literature-based spectral observa- 
tions used to obtain FWHMh/j. Nevertheless, AGNs of the lu- 
minosity considered here usually vary by only ~ 0.13 mag on 
long timescales (e.g., Giveon et al. 1999). The widths of broad 
emission lines in type 1 AGNs typically have an uncertainty of 
~ 10% (e.g., Marziani et al. 2003). Taking all of these factors 
into consideration, we estimate that they introduce a measure- 
ment uncertainty of only ^0.2 dex in Mbh- The largest source 
of uncertainty for the single-epoch masses, however, probably 
comes from the intrinsic scatter of the radius-luminosity rela- 
tion, which is estimated to be ~ 0.4 dex (Bentz et al. 2006). 
According to Peterson et al. (2004), the values of Mbh derived 
from reverberation mapping are accurate to ~ 30%, or 0. 1 dex. 
We note that there is an additional uncertainty on the geometri- 
cal factor (/). For instance, Collin et al. (2006) argued that / 
might be dependent on the accretion rate. We visit this issue in 
§5.2. 

4. THE M BH -ibui RELATION FOR TYPE 1 AGNS 

Figure 5 shows the Mbh _ ^bui relation for our sample of 
type 1 AGNs. Objects with Mbh derived from reverberation 
mapping are encoded separately from those based on single- 
epoch spectra. We assume that the Mbh _ ^bui relation follows 
a log-log relation 



log(M B H/M Q ) = a + /3M s 



bul ■ 



(5) 



The thick hatched line represents the MBH-^bui relation for in- 
active galaxies from Bettoni et al. (2003). Converted to our 
cosmology (see Peng et al. 2006a), the best-fitting relation for 
inactive galaxies is described by a = -2.6 and j3 = -0.5, with a 
scatter of 0.4 dex. 

We estimate a and /3 for the active sample by minimizing \ 2 , 
defined as 



N 

1=1 



(Si 



-(3xd 2 



4 + /3 



2 e 2. 



(6) 
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FIG. 6. — Dependence of the Mbh — Aid relation on Eddington ratio. Blue 
stars and solid line represent the correlation and the best fit for high-Eddington 
ratio objects (Z-bol^Edd^ 0.1). Red circles and dashed line represent low- 
Eddington ratio objects (Lb o i/iEdd<0.1). We fixed the slope to -0.5, as derived 
for the inactive galaxies, whose relation is denoted by the thick hatched line. 
Measurement uncertainties of 0.1 and 0.2 dex are adopted for Mbh estimated 
from reverberation mapping and single-epoch spectra, respectively. 

where y = log(MBH/M Q ), x = Mg^ui, and e,,,- and e X j are the mea- 
surement errors of y andx, respectively (Tremaine et al. 2002). 
This method treats both x and y as independent variables and 
accounts for asymmetric uncertainties for each. 

The estimation of the MBH-^bui relation depends on e Vi , the 
choice of errors for M BH . If we adopt an uncertainty of 0.4 dex 
for the masses based on single-epoch spectra, the x 2 value is 
practically dominated by the reverberation-mapped objects be- 
cause their uncertainties are a factor of 4 smaller, resulting in 
a MBH-^bui relation strongly biased toward the reverberation- 
mapped subsample. For concreteness, we assume that uncer- 
tainties on the single-epoch masses are 0.2 dex, which is a typ- 
ical measurement error (§3.3). As Figure 5 shows, our sample 
of AGNs cluster around the fiducial Mbh _ ^bui relation of in- 
active galaxies with significant scatter. The formal fit for the 
AGNs has a slope of j3 = -0.26 ± 0.05, flatter than for inactive 
galaxies (J3 = -0.5), but because of the limited dynamic range in 
Mbh, we do not regard the AGN fit to be robust. A more mean- 
ingful exercise is to fix the slope of the relation to the value for 
inactive galaxies and then examine the offset and scatter of the 
AGN sample. Fixing to -0.5, the AGN sample has Aa = -0.3 
and an rms scatter of 0.4 dex. 

4. 1 . Dependence on Eddington Ratio 

To understand the physical origin of the intrinsic scatter in 
the M B h _ ^bui relation, we divide the sample into two bins in 
Eddington ratio, at Lboi/^Edd = 0.1. The Eddington luminosity is 
defined as L EdA = 1.26 x 10 38 (M B h/M q ) ergs s" 1 , and the bolo- 
metric luminosity is estimated assuming Lboi = 9ALsioo (Kaspi 
et al. 2000). Figure 6 shows a clear offset between the two sub- 
samples. At a given Mbh, objects with high Eddington ratios 
tend to be hosted by more luminous bulges, or, alternatively, at 
a given bulge luminosity they tend to have less massive BHs. 
In order to quantify the offset between the two subsamples, 
we fix the slope to that of the Mbh _ ^bui relation for inactive 
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FIG. 7. — Similar to Figure 6, except that here we show the dependence on 
galaxy morphology: elliptical galaxies (B/T = 1; blue stars and solid line), 
bulge-dominated systems (0.5 < B/T < 1; green squares and dash-dotted 
line), and disk-dominated systems (B/T < 0.5; red circles and dashed line). 

galaxies. At a fixed M^bui, the offset in Mbh is ~ -0.6 dex; at a 
fixed Mbh, the offset in Mr^uI is ~ 1.2 mag. We note that these 
offsets are much larger than the measurement errors. Perform- 
ing a Kolmogorov-Smirnov test to evaluate how M R ^ in the 
two subsamples is distributed, we find that the null hypothesis 
that the two subsamples are drawn from the same parent popu- 
lation can be rejected with a probability of 97.3%. As discussed 
in §5.2, the segregation between the two subsamples really do 
reflect intrinsic differences in Eddington ratios rather than un- 
certainties in the determination of BH mass. 



4.2. Dependence on Morphological Type 

The availability of robust structural decomposition gives us 
an opportunity to examine possible trends with morphological 
type. Using the measured values of bulge-to-total luminosity 
ratio (B/T; Table 2) and the correlation between morphological 
type and B/T in normal, inactive galaxies (Simien & de Vau- 
couleurs 1986), we divide the sample into three subgroups: 
B/T = 1 (ellipticals), 0.5 < B/T < 1 (bulge-dominated), and 
B/T < 0.5 (disk-dominated). Figure 7 (see also Table 4) shows 
that the zero point of the Mbh -£bui relation, and possibly scat- 
ter, may depend on B/T, although given the limited statistics 
we regard the evidence as tentative. Ellipticals and early-type, 
bulge-dominated galaxies appear virtually indistinguishable, 
but later-type, disk-dominated systems (B/T < 0.5) appear dis- 
tinctly offset to larger Mbh (by ~ 0.4-0.6 dex) at a fixed Mr bui- 
The magnitude of the offset is much larger than possible sys- 
tematic biases in bulge luminosities resulting from uncertainties 
in bulge-to-disk decomposition (~ 0.2 mag; see Fig. 14 in Kim 
et al. 2008). As discussed in §5.1, in many instances our two- 
component fits may not correspond strictly to a bulge+disk de- 
composition but rather to a bulge+tidal feature decomposition. 
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FIG. 8. — Similar to Figure 6, except that here we show the dependence 
on the degree of interaction. Non-interacting or mildly interacting objects 
(ay < 0.3) are denoted by blue stars and the solid line, while strongly inter- 
acting systems (ay > 0.3), with their names labeled, are denoted by red circles 
and the dashed line. Four sources with probable minor companions are marked 
as filled blue stars. The A/bh — ^bul relation for inactive galaxies is shown by 
the thick hatched line. 



4.3. Dependence on Tidal Interaction 

We make use of the quantitative measure of galaxy asymme- 
try, «i, to study the possible effect of tidal interaction. Given 
our small sample size, we simply divide it into two bins accord- 
ing to the value of a\. As shown in Figure 4, a\ rs 0.1 seems 
to provide a useful empirical boundary between objects that are 
disturbed morphologically (a t > 0.1) from those that are not 
(a\ < 0.1). With this threshold for a\, however, the two pop- 
ulations show no obvious segregation in the MfiH-ibui plane. 
But with the boundary set at a higher threshold of a\ = 0.3, Fig- 
ure 8 illustrates that four out of the six objects in our sample 
with the clearest signs of morphological disturbance do prefer- 
entially seem to lie among the most extreme negative outliers 
in the MRH-^bui relation. 

4.4. Dependence on Radio Properties 

The physical drivers responsible for the generation of jets 
and radio emission in AGNs are still unclear. Suggestions have 
included BH mass (e.g., Laor 2000), accretion rate (e.g., Ho 
2002), and host galaxy morphology, which might ultimately be 
linked to BH spin (Sikora et al. 2007). Figure 9 shows that 
a clear separation exists between radio-loud objects and radio- 
quiet objects. Radio-loud sources lie preferentially below the 
Mrh - Lbui relation of inactive galaxies and radio-quiet sources. 

4.5. Dependence on Redshift 

The redshift range of our objects is small (0 < z < 0.35), 
and our sample was not designed to test for evolutionary ef- 
fects. Nevertheless, even by z = 0.36 Woo et al. (2006) and 
Treu et al. (2007) have claimed that AGNs already show evi- 
dence of evolution in the BH-host galaxy scaling relations. Di- 
viding the sources into two bins in redshift (Fig. 10), it ap- 
pears that the two subsamples are offset from each other in 



FIG. 9. — Similar to Figure 6, except that here we show the dependence 
on radio emission: radio-loud sources are marked with blue stars and solid 
line, and radio-quiet sources are marked with red circles and dotted line. The 
Mbh - AjuI relation for inactive galaxies is shown by the thick hatched line. 

the sense that lower-redshift sources have a higher Mrh at a 
given Ms bui- This trend, however, is probably a selection ef- 
fect because low-Eddington ratio, later-type galaxies tend to be 
closer. Indeed, the z < 0.15 subsample has (Lboi/^Edd) =0.11 
and (B/T) = 0.71, to be compared with (Lboi/^Edd) = 0.24 and 
(B/T) = 0.94 for the z > 0.15 subsample. 



5. DISCUSSION 

5.1. Which is the Primary Variable? 

We have assembled a sample of local massive type 1 
AGNs with reliable spectroscopic data and host galaxy 
photometric measurements to investigate the origin of 
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the intrinsic scatter in the correlation between BH mass and 
bulge luminosity. Assuming a geometrical factor of / = 0.75 
for the BLR, we find that the AGNs in our sample lie system- 
atically below the Mbh -£bui relation of inactive galaxies by an 
average offset of Aa « -0.3 dex. Moreover, we have shown 
that the magnitude of the offset correlates with secondary pa- 
rameters connected with the AGN (Eddington ratio and degree 
of radio-loudness) and host galaxy (redshift, bulge-to-disk ra- 
tio, and signs of morphological disturbance). 

Among the several variables that correlate with the offset in 
the Msu-Lbui relation, the only one that appears unphysical is 
that related to redshift. As we noted in §4.5, the apparent de- 
pendence on redshift most likely reflects the selection effect that 
higher redshift sources tend to be biased toward higher Edding- 
ton ratios (e.g., Boyle et al. 2000) and more luminous, earlier 
Hubble types. If we divide the sample into two at z = 0.15, the 
MeH-ibui relation for both the nearby and distant halves con- 
tinue to exhibit the dependence on Lbd/Ledd and B/T that we 
see for the full sample. 

Still, among the rest of the variables that correlate with the 
zero point offset in the Mbh - L\, u i relation, which is the pri- 
mary one? Given that many galaxy and AGN parameters are 
mutually correlated, this is not a trivial question to answer. We 
propose that the primary physical driver is the mass accretion 
rate, as reflected in the Eddington ratio. We argue that the host 
morphology and degree of tidal disturbance directly affect the 
AGN accretion rate, and that the accretion rate, in turn, is linked 
to the radio-loudness parameter. 

Although our 2-D fits indicate that the hosts of many of our 
AGNs have a non-zero disk component apart from a bulge, 
it is important to recognize that, with few exceptions (Fairall 
9, HE 0306-3301, MS 1059.0+7302, PG 2130+99), most of 
the sources in our sample do not have regular, normal disks. 
The vast majority of the sample — by construction when we im- 
posed the Mbh cut — is decisively bulge-dominated. Many of 
the features that we attribute to a "disk," in fact, are simply dif- 
fuse, extended features above and beyond the dominant bulge 
component, which we have parameterized using a single Sersic 
function. There is no a priori reason why the bulge should be 
defined in such a manner, that it cannot have a more complex 
light distribution, especially at large radii. In other instances, 
the extra-bulge component is highly disturbed and almost cer- 
tainly of tidal origin. With few exceptions (Bennert et al. 2008), 
these features have generally never been measured before quan- 
titatively in AGN host galaxies. However, it is entirely debat- 
able whether any of these structures truly belongs to or will 
ever settle into a normal disk component. Instead, we surmise 
that many of the tidal tails and extended, distorted features, in 
fact, should be considered as part of the bulge in formation. 
They are reminiscent of morphological signatures attributed to 
the late, advanced stages of gas-rich mergers (e.g., Barnes & 
Hernquist 1996; Lotz et al. 2008) or possibly even gas-poor 
("dry") mergers (e.g., van Dokkum 2005; Naab et al. 2006). 
Plausible examples of this phenomenon in our sample include 
[HB89] 0316-346 (Fig. 15), HE 1434-1600 (Fig. 19), and 
PG 1012+008 (Fig. 32). In support of this hypothesis, Fig- 
ure lib illustrates that the morphological segregation seen in 
the MfiH-ibui relation (Fig. 7) essentially disappears when the 
bulge luminosity is replaced with the total luminosity of the 
host. The scatter also goes down slightly, from 0.40 dex to 0.36 
dex in the Mbh _ ^host relation. The dependence on Eddington 



ratio, however, remains (Fig. 1 la). 

The above interpretation offers a plausible explanation for 
the apparent link between morphological type and accretion 
rate, which otherwise is somewhat perplexing. Within our sam- 
ple it is the apparently earlier-type, more bulge-dominated sys- 
tems that actually have higher accretion rates. This runs counter 
to the trend normally seen in nearby AGNs (e.g., Heckman 
et al. 2004; Greene & Ho 2007a) and the general tendency 
for present-day early-type galaxies to be more gas-poor than 
late-type galaxies. However, if high-mass, luminous AGNs re- 
sult from the aftermath of gas-rich galaxy-galaxy mergers (e.g., 
Sanders et al. 1988; Hopkins et al. 2006), our results imply 
that it is during the most advanced stages of the merger that ac- 
cretion on the central BH attains its maximum rate. Some of 
the most highly accreting objects (those with large Lboi/^Edd) 
in our sample, which coincide among those with the largest 
offset in the Mbh _ ^bui relation, also happen to be among the 
ones that show the most conspicuous morphological signatures 
of tidal perturbations, as measured by the Fourier parameter 
a\ (see Fig. 8). These include HE 0354-5500 (Lboi/iEdd= 
0.57; a\ = 0.33), PG 1012+008 (L bol /L Edd = 0.32; a x = 0.32), 
PG 1613+658 (L b oi/i E dd= 0.44; a\ = 0.40), and PG 1700+518 
(Lboi/iEdd= 0.71; a\ = 0.39). However, not every object with 
a high Lboi/^Edd has a large value of a\. This may imply that 
not all quasar episodes are triggered by major galaxy interac- 
tions, or that accretion can proceed at a substantial rate even 
after the tidal features have disappeared. To the extent that the 
peak star formation rate in the merger has already subsided dur- 
ing this phase, our scenario offers an additional explanation for 
why type 1 AGNs contain intermediate-age stars (Kauffmann 
et al. 2003) but generally not much concurrent star formation 
(Ho 2005; Kim et al. 2006). 

Within our sample, four AGNs have compact sources — 
plausibly small accreted companions — projected close to the 
primary host galaxy. These may be examples of minor mergers. 
None has a large value of a\ . Two of the four lie exactly on the 
MBH-^bui relation of inactive galaxies (Fig. 8) and have rela- 
tively low Eddington ratios (L bo i/L E dd = 0.05 for PG 1426+015 
and L hol /L EM = 0.06 for PKS 2135-14). The other two lie offset 
below the relation, but only PG 1302-102 has a large Lboi/^Edd 
(0.4); PHL 1093 has L hol /L EM = 0.03. Thus, within our lim- 
ited statistics, we have no evidence that minor mergers play a 
significant role in elevating the accretion rate in AGNs. 

Within this backdrop, we can offer a tentative explanation for 
the zero point difference between radio-loud and radio-quiet 
sources (Fig. 9), one that ultimately links the generation of 
powerful radio jets to the BH accretion rate and/or host galaxy 
morphology. But first we should clarify some terminology. 
There are two popular definitions of "radio-loud" AGNs in the 
literature, and it is important not to confuse them. One common 
usage of this term refers to sources that are classified solely by 
their radio-to-optical flux ratio (R) as defined by Kellermann et 
al. (198913. regardless of their radio power. On this basis, the 
vast majority of AGNs in the local Universe (Ho 2008), most 
with extremely low luminosities, qualify as being radio-loud, 
with the degree of radio-loudness increasing with decreasing 
Lboi/^Edd (Ho 2002; Terashima & Wilson 2003; Greene et al. 
2006). The radio emission in most of these low-power sources 
is largely dominated by a compact core, and any jet-like fea- 
tures are confined to sub-galactic scales. The host galaxies en- 
compass all morphological types, including spiral galaxies (Ho 



Terashima & Wilson (2003) advocated a closely related radio-loudness parameter Rx based on the radio-to-X-ray flux ratio. 
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between Mbh and bulge luminosity for inactive galaxies is denoted by the thick hatched line. 



& Peng 2001). The Milky Way's Sgr A* is a familiar example. 
The second commonly used definition of radio-loudness is less 
clear-cut, but it involves some combination of relative (R ^> 10) 
and absolute (P m d ;> 10 23-24 W Hz" 1 ) measures of radio power. 
When detected, the radio jets have super-galactic dimensions 
and are highly collimated. The radio-loud sources in this study, 
with median R = 820 and P^ = 7 x 10 25 W Hz" 1 , belong to 
this second category. Strong radio sources of this variety in- 
variably reside in early-type galaxies (e.g., McLure et al. 1999) 
and are associated with high accretion rates (e.g., Maccarone et 
al. 2003; Kording et al. 2006). Nevertheless, only a minority 
of highly accreting AGNs in massive, early-type host galaxies 
are radio-loud. There is no universally accepted explanation for 
this longstanding quandary. One possibility is that a necessary 
ingredient for the generation of powerful, collimated jets is the 
existence of a BH with a large spin (e.g., Sikora et al. 2007), 
which is more easily attained in merger-driven accretion events 
during the formation of giant elliptical galaxies than in disk or 
spiral galaxies (Volonteri et al. 2007). 

5.2. Physical Interpretation 

The principal conclusions of our study are that the zero point 
of the Mbh -ibui relation for AGNs is offset from that of inac- 
tive galaxies and that the magnitude of the offset correlates with 
several physical properties of the AGN and host ultimately con- 
nected to the accretion rate. These effects account for the bulk 
of the observed intrinsic scatter in the AGN MBH~£bui relation. 

Interestingly, other studies of nearby AGNs using different 
probes of the host galaxy have independently arrived at very 
similar conclusions. Onken et al. (2004; see also Nelson et al. 
2004), analyzing a small sample of type 1 AGNs with available 
stellar velocity dispersions and BH masses determined through 
reverberation mapping, concluded that the active systems lie 
offset below the Mbh _ c* relation of inactive galaxies by ~ 0.2 
dex if one assumes a spherical BLR. Greene & Ho (2006b) ob- 



tained an almost identical result using a much larger sample of 
type 1 AGNs with measured a* and Mbh estimated from single- 
epoch spectra. These conclusions have been reaffirmed by Shen 
et al. (2008) in their analysis of composite AGN spectra derived 
from the Sloan Digital Sky Survey. A different approach was 
taken by Ho et al. (2008) through 21 cm H I observations. Us- 
ing the maximum rotational velocity and total dynamical mass 
of the galaxy as new variables to represent the gravitational po- 
tential of the host, they find that both quantities strongly corre- 
late with Mbh- In qualitative agreement with the results from 
this study, Ho et al. find that the zero point in the scaling rela- 
tions depends primarily on the accretion rate in the sense that, at 
a given galaxy rotation velocity and especially dynamical mass, 
AGNs with higher Eddington ratios have systematically lower 
Mbh than those with lower Eddington ratios. 

These trends can be interpreted in one of two ways. In the 
first instance, we note that AGNs, by selection, have actively 
growing BHs. We can envision that, at a fixed bulge poten- 
tial (velocity dispersion, luminosity, mass), active galaxies have 
less massive BHs than inactive galaxies if the bulk of the star 
formation precedes, and is not well synchronized with, a ma- 
jor accretion event. This particular time sequence, which is re- 
quired in order to imprint a net negative offset in the Mbh versus 
host galaxy relations, seems to be supported by the prevalence 
of post-starburst signatures (Canalizo & Stockton 2001; Kauff- 
mann et al. 2003) as well as the low levels of ongoing star 
formation (Ho 2005; Kim et al. 2006) found in type 1 AGNs. If 
the AGN phase lasts, say, for ~ 10 8 yrs, which is near the upper 
end of the currently estimated lifetimes (Martini 2004), then a 
10 8 Mq BH would increase its mass by a factor of 2 (~ 0.3 dex) 
if it is radiating at Lboi/^Edd = 0.5 with a radiative efficiency of 
0.1. This example is merely illustrative. In reality, the AGNs in 
our sample span a wider range of Lboi/^Edd and the AGN life- 
times may be shorter. Nevertheless, on average we expect lumi- 
nous AGNs to lie systematically below inactive galaxies on the 
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Mbh~L\,uI relation. Moreover, our study, as do those of Ho et 
al. (2008) and Shen et al. (2008), further shows that the magni- 
tude of the zero point offset depends on the Eddington ratio: the 
BH masses of high-Eddington ratio AGNs have more catching 
up to do than the BH masses of lower-Eddington ratio AGNs. 
Such a systematic trend can only come about if the accretion 
rate directly relates to the evolutionary phase of the accretion 
event. This seems plausible, in light of the apparent association 
between accretion rate and the degree of tidal perturbation. 

Instead of the BH being undermassive, perhaps it is actually 
the bulge that is overluminous, by AM«,bui ~ 0.5 - 0.6 mag. 
First, we dismiss the possibility that the luminosity enhance- 
ment could be due to contamination from nebular emission 
from the narrow-line region. Although the spatial extent of 
the narrow-line region in quasars can reach several kpc (Ben- 
nert et al. 2002), substantially overlapping with the bulge, 
the typical [O III] luminosities in our sample ((L[om]) = 10 42 
ergs s" 1 ) contribute less than 3% to the luminosity of the bulge 
((Mftbui) = -22.75 mag). Given the evidence outlined in §5.1 
that the most extreme outliers in the Mbh - ^bui seem to have 
undergone a recent merger or tidal interaction, a more likely 
possibility is that the bulge luminosity may be moderately en- 
hanced by the latest episode of central star formation. Indeed, 
for a small sample of reverberation-mapped Seyfert 1 galaxies 
with stellar velocity dispersion and bulge luminosity measure- 
ments, Nelson et al. (2004) have shown that these objects are 
somewhat brighter (~ 0.4 mag) than inactive galaxies at a given 
velocity dispersion. They attributed this offset in the Faber- 
Jackson (1976) relation to younger stellar populations in AGNs. 
Boosting the /?-band luminosity by 0.5 mag requires ~ 15% of 
the stellar mass to come from a 1 Gyr population with solar 
metallicity (Nelson et al. 2004). While this offers a plausible 
explanation for the offset in the Mbh-£|ju1 relation seen in our 
sample, it cannot account for the fact that the most recent and 
largest samples of nearby type 1 AGNs, statistically at least, 
show negative offsets (with respect to inactive galaxies) when 
Mbh is compared to all bulge or host galaxy parameters (stellar 
velocity dispersion, bulge luminosity, total galaxy dynamical 
mass). The direction of the offset is the same (at a given galaxy 
parameter, Mbh is lower), and the magnitude of the offset is 
also similar (~ 0.2-0.3 dex). Although this clearly needs to be 
verified with a large sample that has reliable measurements of 
both velocity dispersion and bulge luminosities, the most recent 
studies suggest, contrary to Nelson et al. (2004), that local type 
1 AGNs actually do not depart from the standard Faber-Jackson 
(1976) relation. Furthermore, for a large sample of sources with 
measurements of rotational velocity and total galaxy luminos- 
ity, Ho et al. (2008) show that type 1 AGNs show no obvi- 
ous deviations from the Tully-Fisher (1977) relation of inactive 
galaxies. In light of these considerations, we favor the view that 
the negative offset in the MBH-^bui relation represents a deficit 
in Mbh rather than an excess in Lbui- 

Alternatively, we can assert that both active and inactive 
galaxies intrinsically should obey the same BH mass versus 
host scaling relations. From this standpoint, the zero point off- 
set between AGNs and inactive galaxies, as well as the varia- 
tions of the offset with Lboi/^Edd, can be viewed as a system- 
atic underestimate of the true value of Mbh- Recall that our 
BH masses are based on a virial product assuming a spherical 
distribution of BLR clouds with isotropic velocities, for which 
the geometric factor is / = 0.75. If, for example, the BLR (or 
at least the portion of it that predominantly emits the Balmer 



lines) has a flattened, disk-like geometry with kinematics dom- 
inated by rotation, and on the scale of the BLR type 1 sources 
happen to be preferentially more face-on to our line of sight, 
then we systematically underestimate the deprojected rotation 
velocity and hence Mbh- Wu & Han (2001) invoked this line of 
reasoning to interpret the offset of type 1 AGNs on the Mbh-o* 
relation and concluded that on average their BLRs are inclined 
by (/) w 36°. The same argument can, in principle, be applied 
to the observed offset in the MBH-^bui relation. While the ef- 
fect of inclination probably enters at some level (see also Collin 
et al. 2006), it cannot account for the fact that the magnitude of 
the offset depends on Lboi/^Edd- The latter is not a trivial conse- 
quence of the mass being underestimated because sources with 
high Lboi/^Edd truly do exhibit characteristically distinct X-ray, 
optical, and radio properties (Boiler et al. 1996; Boroson 2002; 
Greene et al. 2006). Whatever the physical origin of the offset 
(Onken et al. 2004; Collin et al. 2006; Marconi et al. 2008), 
we can empirically adjust the normalization factor of the virial 
product by forcing the AGN sample to agree with the fiducial 
reference of inactive galaxies. To remove the zero point offset 
of Aa w -0.2 to -0.3 dex, then, the normalization should be in- 
creased by a factor of ~ 1.6-2, from/ = 0.75 to/« 1.2-1.5. 
For the most extreme offsets of Aa s=s -0.6, / increases to ~ 3. 

Without additional information, unfortunately, the above two 
alternative explanations — an undermassive Mbh versus an un- 
derestimated / factor — are degenerate. It is easy to imagine 
that both effects must operate jointly. On the one hand, the 
BHs in AGNs are, after all, gaining mass. On the other hand, as 
discussed in Collin et al. (2006), there are multiple reasons to 
believe that the BLR has a nonspherical geometry and that the 
Eddington ratio may influence its structure and dynamics. The 
only way to resolve this degeneracy is to obtain independent 
estimates of Mbh for AGNs that do not rely on the BLR virial 
technique. To date, efforts to apply resolved stellar dynamical 
techniques to reverberation-mapped AGNs have yielded very 
rough estimates of Mbh for only a couple of sources (NGC 
3227: Davies et al. 2006; Hicks & Malkan 2008; NGC 4151: 
Onken et al. 2007), and thus attempting to cross-calibrate the 
two techniques is still far too premature. BH mass estimators 
based on X-ray variability seem more promising. Gierlihski et 
al. (2008; see also Hayashida et al. 1998) find that, for accret- 
ing BHs in their hard spectral state, the amplitude of their high- 
frequency X-ray variability scales inversely with Mbh over a 
very wide range of masses. For a small subset of nearby 
type 1 AGNs, the X-ray-derived masses show rough agreement 
with Mbh obtained through reverberation mapping assuming 
/ s=a 1 .2. A very similar conclusion was reached by Nikolajuk et 
al. (2006). From comparison of Mbh for reverberation-mapped 
sources with masses obtained using the X-ray excess variance 
method, these authors estimate /= 1.06 ±0.26, which, inter- 
estingly, lies in between the values of the / factor for the two 
extreme alternatives discussed above. 

Assuming that the X-ray-derived normalization factor truly 
does represent the correct normalization factor for the virial 
masses, then the inferred growth rates for Mbh in luminous 
AGNs are much more modest, from typically as little as 10%- 
40% (0.05-0.15 dex) to at most 280% (0.45 dex). 

Still, we note that the tendency for BHs in local AGNs to be 
less massive than the BHs in inactive galaxies of similar type 
runs counter to the trend observed at higher redshift. Already by 
z = 0.36, type 1 AGNs begin to depart from the local Mbh - c* 
and M BH -ibui relations (Woo et al. 2006; Treu et al. 2007), 
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but in the opposite direction as that seen at lower redshifts: at 
a given a+ or Lt, u i, AGNs are offset compared to local inactive 
systems by AlogM BH ~ +0.5 dex. This trend has now been ex- 
tended by Woo et al. (2008) out to z = 0.57 using stellar velocity 
dispersion measurements. At even higher redshifts, direct a* 
measurements are no longer feasible, but surrogate estimates of 
cr* using narrow emission lines (Salviander et al. 2007) as well 
as probes of the host galaxy using imaging (Peng et al. 2006a, 
b) and CO emission lines (Shields et al. 2006; Ho 2007) sup- 
port the notion that the growth of the BHs in AGNs have been 
decoupled from, and outpaced, the underlying host. 

6. SUMMARY 

We performed two-dimensional structural decomposition of 
a sample of 45 nearby (z <^ 0.35) type 1 AGNs with available 
archival optical HST images and published spectroscopic data. 
We calculated virial BH masses assuming a spherical BLR with 
isotropic velocities. Using a new version of the versatile code 
GALFIT, we derived detailed fits to the structural components 
of the host galaxies, yielding not only robust measurements of 
the bulge luminosities with realistic error bars but also, for the 
first time, quantitative estimates of nonaxisymmetric features 
such as extended disks and tidal arms. 

Our principal aim is to understand the origin of the intrinsic 
scatter in the MsH-Lbui of active galaxies over a restricted range 
of BH masses (M BH ~ 10 85±a5 M Q ). While AGNs closely fol- 
low the Mbh -L\, a i relation of inactive galaxies, we find that the 
intrinsic scatter is substantial (0.40 dex) and that the zero point 
of the relation is shifted by an average of AlogMBH ~ — 0.3 
dex. The magnitude of the zero point offset in the Mbh _ ^bui 



relation depends on properties of the AGN (Eddington ratio and 
radio-loudness parameter) and the host galaxy (morphological 
type and degree of tidal perturbation). We argue that the prin- 
cipal physical parameter responsible for the variation in zero 
point is the BH accretion rate, as reflected in the Eddington ra- 
tio. We suggest that galaxy mergers and tidal interactions play 
a substantial role in boosting the accretion rate in this sample 
of AGNs. A significant fraction of the zero point offset in the 
Mbh _ ^bui relation can be explained if the virial BH masses 
have been underestimated, as indicated from comparison with 
independent masses derived from X-ray variability techniques. 
After accounting for this change in the normalization of the 
virial BH mass scale, we estimate that BHs during the AGN 
phase experience a modest growth of ~10%-40% in mass. 
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APPENDIX 

NOTES ON INDIVIDUAL OBJECTS 
Comments on the fitting results for individual objects are given here. 



3C 59 (Fig. 12) — The host can be fit with classical bulge represented by a de Vaucouleurs (n = 4) profile. 

E 1821+643 (Fig. 13) — The best fit shows that the bulge is slightly disturbed (ci\ = 0.17). There is no evidence of a disk. 

Fairall 9 (Fig. 14) — The host requires a bulge and a disk. A ring-like structure, which we do not model, is also seen in the 
residuals. The bulge appears be quite compact, with an effective radius of 1.6 pixels, although it is not well decomposed from the 
nucleus. Thus, the bulge luminosity might be highly uncertain. 

[HB89] 0316-346 (Fig. 15) — Highly disturbed object with prominent tidal features. We fit the host with an n = 4 bulge and two 
disk components with Fourier modes. One of the disk components is not centered on the nucleus. 

[HB89] 2201+315 (Fig. 16) — This object was observed in two different filters (F555W and F702W), but both images are 
saturated in the core. The bulge luminosities derived from the two different images (corrected to the R band) differ by 0.7 mag. We 
fit the host with a single n = 4 bulge in both images. The residual image from the long (2700 s) exposure shows possible signs of an 
extended disk. 

HE 0306-3301 (Fig. 17) — The best-fit result shows an elongated bulge component (b/a s=s 0.45) and a disk with a spiral arm. 
The signal-to-noise ratio (S/N) of the observed PSF star is lower than that of the science image. 

HE 0354-5500 (Fig. 18) — This appears to be a merging system. We fit the host galaxy with a bulge (n = 4) and an off-centered 
tidal-like feature with a high-amplitude Fourier mode {a\ = 0.33). The S/N of the observed PSF star is lower than that of the science 
image. 

HE 1434-1600 (Fig. 19) — We fit the host with only a single bulge component (m\, u \ w 17.5 mag), although the residual image 
shows evidence for arcs and ripples, which may be evidence of a recent collision (Letawe et al. 2004). Alternatively, if the host is fit 
with two components, the bulge magnitude becomes 18.6 mag. The S/N of the observed PSF star is lower than that of the science 
image. 

MC 1635+119 (Fig. 20) — We fit the host with a single bulge component (n = 4). 

MRK 1048 (Fig. 21) — This object has a close companion and an extremely large but faint, off-center tidal feature. However, 
neither of these has a large effect on the fitting result. We fit the host with a classical bulge (n = 4) and a tidal tail. 

MS 0244.8+1928 (Fig. 22) — The host is reasonably well represented by an n = 4 bulge, although the residuals indicate that 
there might be an additional faint, outer envelope. 

MS 0754.6+3928 (Fig. 23) — The host is reasonably well represented by an n = 4 bulge, although the residuals indicate that 
there might be an additional faint, outer envelope. 
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MS 1059.0+7302 (Fig. 24) — The best-fitting model requires a bulge and a disk. There is a ring-like structure in the residual 
image. 

MS 1545.3+0305 (Fig. 25) — The best fit requires a bulge (n = 4) and a disturbed exponential disk. A ring-like structure is 
prominent in the original and residual images. 

OX 169 (Fig. 26) — There are significant residuals after subtracting the best-fit observed PSF, which might be due to PSF 
mismatch. The fit with the Tiny Tim PSF results in a brighter AGN by 0. 1 mag. The host contains a bulge (n = 3.9) and a prominent, 
extended, asymmetric structure that is fit with a Sersic component with a small n (~ 0.2). 

PG 0052+251 (Fig. 27) — It appears to have a tidally disturbed spiral arm, which may be related to the two small galaxies in its 
vicinity. We fit this object with a classical bulge (n = 4) and a truncated disk with a small n (~ 0.2). 

PG 0804+761 (Fig. 28) — There appears to be a faint central feature that resembles a bar or highly inclined disk-like structure. 
The bulge luminosity, however, is hardly affected by the bar component. 

PG 0923+201 (Fig. 29) — The host is fit with a single bulge component (n = 4). There are three nearby companions that were 
fit simultaneously. 

PG 0953+414 (Fig. 30) — The host is fit with a single bulge component (n = 4), which appears slightly disturbed based on the 
amplitude of the Fourier mode (ai = 0.12). 

PG 1004+130 (Fig. 31) — The host is fit with a single bulge component (n = 4). The residual image shows significant PSF 
mismatch, which might affect the fitting result. 

PG 1012+008 (Fig. 32) — The image of the host shows clear evidence of interaction with a spiral galaxy and a smaller, compact 
neighbor. The primary host can be fit with a single bulge component (n = 4), but the residual image shows significant structure. 

PG 1116+215 (Fig. 33) — The residual image indicates that the central core is slightly affected by PSF mismatch. The host 
galaxy, however, is well-represented by a single-component bulge with n fixed to 4. If we allow n to be free, the best-fitting value of 
n = 1 .84 yields a bulge luminosity that is 0.3 mag fainter. 

PG 1202+281 (Fig. 34) — This source has a number of nearby galaxies, the brightest and nearest of which we fit simultaneously. 
The primary host is well described by a single-component bulge (« = 4). 

PG 1211+143 (Fig. 35) — Like PG 0804+761, there appears to be a faint central feature that resembles a bar or highly inclined 
disk-like structure. It is unclear if this is an artifact due to PSF mismatch. Depending on whether this extra component is included, 
the bulge luminosity ranges between 16.7 and 17.2 mag, with a best-fit value of 16.9 mag. 

PG 1226+023 (Fig. 36) — Because of the strong bleeding regions from the saturated core, we did not perform nonparametric 
aperture photometry. The host is reasonably well fit with a single bulge component (n = 4). The residual image shows the well-known 
jet of this source (3C 273), as well as some diffuse, extended emission. The box-like imprint in the residual image results from the 
PSF image being smaller than the science image. 

PG 1302-102 (Fig. 37) — This object was observed in two different filters (F606W and F702W), but both images are saturated 
in the core. The image contains two compact sources superposed on the main host, which we fit simultaneously. The host can be fit 
with a single bulge component, although significant structure on large scales remain in the residual image. The fits from both filters 
are in good agreement to within the uncertainty. 

PG 1307+085 (Fig. 38) — Although the core of the image is saturated, the host is well-described by a single bulge component 
(n = 4). 

PG 1309+355 (Fig. 39) — The core of the image is saturated. We fit the host with a single bulge component (n = 4). The residual 
image shows evidence of spiral-like substructure, but we did not attempt to model it. 

PG 1351+640 (Fig. 40) — This is an almost face-on system with spiral arms. The fit is done with a classical bulge (n = 3.8) and 
an exponential disk. 

PG 1411+442 (Fig. 41) — This is an extremely disturbed object that appears to have a nearby companion. We fit this object with 
a classical bulge (n = 4) and a spiral disk with Fourier modes. 

PG 1416-129 (Fig. 42) — The host is well-described by a single bulge component (n = 4). 

PG 1426+015 (Fig. 43) — This object has a tidal tail and a small companion. The best-fitting model for the host consists of a 
pseudobulge (n = 2. 1) and a disk with Fourier modes. Modeling the host with only a single component yields a much brighter bulge 
(14.2 vs. 16.08 mag), but the residuals of the fit are significantly worse than those of the two-component fit. 

PG 1444+407 (Fig. 44) — The host is fit with a bulge (n = 4) and a somewhat disturbed disk component (n = 0.31, a\ = 0.10). 
The residuals suggest that a ring-like component might be present. 

PG 1613+658 (Fig. 45) — This is a highly disturbed object in a merging system. The fit is ambiguous. The best-fitting model 
for the host consists of a classical bulge (n = 4) and a disturbed disk (n = 1) with Fourier modes. The bulge is ~0.7 mag fainter than 
the best-fitting case if the fit is done with a single bulge component (n = 4), but the residuals of the fit are significantly worse than 
those of the two-component fit. 

PG 1617+175 (Fig. 46) — Although the best-fitting model for the host contains a bulge and a disk, a single-component model 
also works reasonably well. In the two-component fits, the bulge luminosities range from 17.7 mag (n = 4) to 17.9 mag (n = 1). 

PG 1700+518 (Fig. 47) — The host is well-represented by a single bulge component (n = 4) plus a tidal tail, which is fit with 
Fourier modes. 

PG 2130+099 (Fig. 48) — This object is very similar to PG 0052+251; both have ring-like spiral disk. Since the central few 
pixels are saturated, the fit is slightly uncertain. The model for the host consists of a pseudobulge (n = 0.44) and a disk (n = 0.33) 
with Fourier modes. 

PHL 909 (Fig. 49) — The host is fit with a classical bulge (n = 4), but the residual image shows a faint central feature that 
resembles a tiny bar or highly inclined disk-like structure. The host galaxy is slightly disturbed. 
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PHL 1093 (Fig. 50) — The host is well fit with a single bulge component (n = 4). Three small nearby companions are included 
in the fit simultaneously. 

PKS 0736+01 (Fig. 51) — The host is well fit with a single bulge component (n = 4). There are several faint blobs nearby, but it 
is unclear whether these are associated with the primary host. 

PKS 1020-103 (Fig. 52) — A single bulge component (n = 4) is enough to describe the host. 

PKS 1217+02 (Fig. 53) — The host is well fit with a single bulge component (n = 4). 

PKS 2135-14 (Fig. 54) — We fit the host with a single bulge component (n = 4). Two nearby objects are included simultaneously 
in the fit. 

PKS 2349-01 (Fig. 55) — There are two curved tidal tails. The best-fit model for the host consists of a classical bulge (n = 4) 
and a highly distorted disk with n « 0.88. 

PKS 2355-082 (Fig. 56) — The host is well fit with a single bulge component (n = 4). 
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FIG. 12. — GALFIT decomposition for 3C 59. (a) Azimuthally averaged profile, showing the original data (open circles), the best fit (solid line), and the 
sub-components (PSF and bulge; dashed lines). The residuals are plotted on the bottom. We present the 2-D image of the original data (b), the best-fit model for the 
host (the AGN is excluded to better highlight the host), with the amplitude of the first Fourier mode (a\) labeled (c), and the residuals (d). The units of the images 
are in arcseconds. All images are on an asinh stretch. 
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FIG. 13. — GALFIT decomposition for E 1821+643; symbols and conventions as in Figure 12. 
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FIG. 14. — GALFIT decomposition for Fairall 9; symbols and conventions as in Figure 12. 



26 



KIM et al. 



-, ||£-:-.w|..'i|,yi4J 1 |".' 




■ ■ 



■ m W- M ■ ■ ■ ■ M 



■I 




m ■ 




k]l«.-*Juj" 




I 



.■ I I 



I T- 



FIG. 15. — GALFIT decomposition for [HB89] 0316—346; symbols and conventions as in Figure 12. 
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FIG. 16.— GALFIT decomposition for [HB89] 2201+315 (PC/F555W); symbols and conventions as in Figure 12. 
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FIG. 17. — GALFIT decomposition for HE 0306—3301; symbols and conventions as in Figure 12. 
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FIG. 18. — GALFIT decomposition for HE 0354-5500; symbols and conventions as in Figure 12. 
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FIG. 19. — GALFIT decomposition for HE 1434—1600; symbols and conventions as in Figure 12. 
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FIG. 20. — GALFIT decomposition for MC 1635+1 19; symbols and conventions as in Figure 12. 
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FIG. 2 1 . — GALFIT decomposition for MRK 1048; symbols and conventions as in Figure 12. 
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FIG. 22. — GALFIT decomposition for MS 0244.8+1928; symbols and conventions as in Figure 12. 
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FIG. 23. — GALFIT decomposition for MS 0754.6+3928; symbols and conventions as in Figure 12. 
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FIG. 24. — GALFIT decomposition for MS 1059.0+7302; symbols and conventions as in Figure 12. 
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FIG. 25. — GALFIT decomposition for MS 1545.3+0305; symbols and conventions as in Figure 12. 



Black Hole Mass vs. Bulge Luminosity Relation 



37 






■ 


4dlPn*hml 




- 






- 


- 




• 


r H a . 








■ 


- 













II k I 






■ 1 I* 



FIG. 26. — GALFIT decomposition for OX 169; symbols and conventions as in Figure 12. 
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FIG. 27. — GALFIT decomposition for PG 0052+251; symbols and conventions as in Figure 12. 
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FIG. 28. — GALFIT decomposition for PG 0804+761; symbols and conventions as in Figure 12. 
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FIG. 29. — GALFIT decomposition for PG 0923+201; symbols and conventions as in Figure 12. 
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FIG. 30. — GALFIT decomposition for PG 0953+414; symbols and conventions as in Figure 12. 
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FIG. 3 1 . — GALFIT decomposition for PG 1004+130; symbols and conventions as in Figure 12. 
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FIG. 32. — GALFIT decomposition for PG 1012+008; symbols and conventions as in Figure 12. 
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FIG. 33. — GALFIT decomposition for PG 1 1 16+215; symbols and conventions as in Figure 12. 
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FIG. 34. — GALFIT decomposition for PG 1202+281; symbols and conventions as in Figure 12. 
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FIG. 35 . — GALFIT decomposition for PG 121 1+143; symbols and conventions as in Figure 12. 
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FIG. 36. — GALFIT decomposition for PG 1226+023; symbols and conventions as in Figure 12. 
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FIG. 37. — GALFIT decomposition for PG 1302—102 (PC/F702W); symbols and conventions as in Figure 12. 
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FIG. 38. — GALFIT decomposition for PG 1307+085; symbols and conventions as in Figure 12. 
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FIG. 39. — GALFIT decomposition for PG 1309+355; symbols and conventions as in Figure 12. 
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FIG. 40. — GALFIT decomposition for PG 1351+640; symbols and conventions as in Figure 12. 
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FIG. 41 . — GALFIT decomposition for PG 141 1+442; symbols and conventions as in Figure 12. 



Black Hole Mass vs. Bulge Luminosity Relation 



53 






— T1 

rfir.-HihHl 




i 


■ 
■ 


V- 




■ 




i ■• 

■ 

■ 

■ 

« 





■ 



I 



FIG. 42. — GALFIT decomposition for PG 1416—129; symbols and conventions as in Figure 12. 
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FIG. 43. — GALFIT decomposition for PG 1426+015; symbols and conventions as in Figure 12. 
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FIG. 44. — GALFIT decomposition for PG 1444+407; symbols and conventions as in Figure 12. 
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FIG. 45. — GALFIT decomposition for PG 1613+658; symbols and conventions as in Figure 12. 
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FIG. 46. — GALFIT decomposition for PG 1617+175; symbols and conventions as in Figure 12. 
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FIG. 47. — GALFIT decomposition for PG 1700+518; symbols and conventions as in Figure 12. 
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FIG. 48. — Example of GALFIT decomposition for PG 2130+099; symbols and conventions as in Figure 12. 
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FIG. 49. — GALFIT decomposition for PHL 909; symbols and conventions as in Figure 12. 
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FIG. 50. — GALFIT decomposition for PHL 1093; symbols and conventions as in Figure 12. 
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FIG. 5 1 . — GALFIT decomposition for PKS 0736+01 ; symbols and conventions as in Figure 12. 
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FIG. 52. — GALFIT decomposition for PKS 1020-103; symbols and conventions as in Figure 12. 
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FIG. 53. — GALFIT decomposition for PKS 1217+02; symbols and conventions as in Figure 12. 
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FIG. 54. — GALFIT decomposition for PKS 2135-14; symbols and conventions as in Figure 12. 
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FIG. 55. — GALFIT decomposition for PKS 2349-01; symbols and conventions as in Figure 12. 
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FIG. 56. — GALFIT decomposition for PKS 2355-082; symbols and conventions as in Figure 12. 



